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[1] Realitatsnaher Mathematikunterricht mit DERIVE, H.-J. Henn, (Diimmlerbuch 4565), Diimmler Ver- 
lag Bonn, ISBN 3-427-4565 1-X. 


Arbeitsblatter zur Mathematisierung realer 
Probleme: Vom Regenbogen iibers Auto und 
Autobahntrassen zu wirtschaftlichen Proble- 
men. 


Working sheets for modelling real life problems: 
the rainbow, some problems In connection with 
cars, how to plot a road, economy and money 
problems. 


[2] Empfehlungen fiir den Mathematikunterricht an Gymnasien, Niedersachsisches Kultusministerium, 
Referat Presse und Offentlichkeitsarbeit, Schilfgraben 12, 30159 Hannover 

Leitlinien zur Themenauswahl und Unterrichtsgestaltung, Zu Veranderungen des Unterrichts, 5 exemplari- 
sche Unterrichtsbeispiele (durchgangig fiir den Einsatz des Tl-92). 


Keith Eames produced a booklet on DERIVE for Windows on Algebra which he uses for his year 12 
students. If you are interested then please contact me or him. You can find his email address in this 
DNL. 

Wilhelm Weiskirch and Heiko Knechtel are asking for partners in the frame work of an EU- 
Comenius project. They both are teachers at a Gymnasium in Niedersachsen, Germany, and they 
want to exchange experiences and to develop supporting materials for upper secondary level teaching 
(special emphasis will be given on the Tl-92). If you are interested then please contact 
w.weiskirch@t-online.de or HKnechtel@aol.com. 


Interesting WEB sites: 

All websites are valid in 2012! 

http://www.maa.org/mathland/mathtrek 2 2 98.html (Ivar Peterson’s MathTrek) 
http://sprott.physics.wisc.edu/pickover/ Both web sites are recommended by Jerry Glynn 

http://metric.ma.ic.ac.uk/came/ The home page of the IC-CAME (International Council for 

Computer Algebra in Mathematics Education) 

http://sunsite.utk.edu/math archives/.http/contests/ Mathematics Contests, Competitions, 

and Problem Sets 

http://ticalc.org A very informative site for Tl-92 users 

http://www.t3ww.org/ The T^ site (Teachers Teaching with Technology) 


Internationales Symposium zur Didaktik der Mathematik 

Thema: Mathematische Bildung und neue Technologien 

Universitat Klagenfurt, 28. Sept. - 2.0kt. 1998 

Information, Application and Submission of papers: 

http://www.uni-klu.ac.at/groups/math/didaktik/symp98/symp98.htm 

or: Universitat Klagenfurt,Inst.f.Mathematik, Abteilung f. Didaktik der Mathematik, 

Universitatsstr. 65, A-9020 Klagenfurt, Tel.: ++43 (0)463 2700 429 
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Liebe DERIVE und TI-Freunde, 

Wieder ist ein Newsletter fast fertig. Es fehlt nur noeh der 
"Letter". Trotz der Arbeit, die es maeht alles in eine ge- 
meinsame Form zu gieBen ist es immer wieder eine Freude, 
eine Ausgabe zusammenzustellen da Hire vielen Beitrage 
das weite Spektrum der DERIVE und TI-92 Anwendungen 
zeigen. Es ist immer eine Herausforderung, die Beitrage fiir 
einen Newsletter zu sammeln, so dass alle in 46 Seiten 
Platz findet und aueh ein wenig zusammenpaBt. 

Wolfgang. Propper hat eine Kolumne fur „Anfanger“ ver- 
sproehen , muB sieh aber fur diesen DNL wegen allzuvieler 
Verpfliehtungen (Vortrage, Kurse) entsehuldigen. leh glau- 
be aber, dass Peter Mities Uberlegungen zum Programmie- 
ren mit DERIVE gerade aueh fur nieht so geiibte DERIVE- 
Anwender von Nutzen sein konnen. Peter hat zwar im Xitel 
"DERIVE 3.x", sein 
Aufsatz hat aber 
naeh meiner 

Meinung aueh in 
Zeiten von DfW 
4.09 und DfD 4.x 
niehts von Ak- 
tualitat eingebiiBt. 

leh maehe Sie geme 
aueh hier auf die 
Suehe von Willi 
Weiskireh und 
Heiko Kneehtel 
naeh Partnem fur 
ein EU-COME- 
NIUS Projekt 
aufmerksam. leh darf Ihnen aber aueh gleiehzeitig berieh- 
ten, dass wir Ende Februar ein Projekt "Modeme Teehno- 
logien im Mathematikunterrieht" bei der EU eingereieht 
haben. Dem Padagogisehen Institut Niederosterreieh als 
Koordinierende Einriehtung haben sieh das ZKL Munster 
(Deutsehland), das CTM an der Uni Plymouth (England), 
die Uni Goteborg (Sehweden), das Freudenthal Institut 
(Niederlande) und die Katholisehe Uni Leuven (Belgien) 
als Partner angesehlossen. leh werde Sie iiber die Entwiek- 
lung weiter informieren - wenn das Projekt bewilligt wird. 

AbsehlieBend moehte ieh Sie auf die beigelegten Program- 
me fur die Gettysburg Conferenee und fur die Miinsteraner 
Pfingsttagung hinweisen. 

Mit vielen GriiBen 


Dear DERIVE and Tl-friends 

And so another Newsletter is nearly ready to go to 
print. Despite the work involved in finding a common 
shape for ail the various contributions each time, I 
have again enjoyed compiling an issue because your 
many contributions are show such a wide field of 
DERIVE and TI-92 applications. It is always a chal- 
lenge to gather articles in a way so that they will fill 46 
pages and that they more or less fittogether. 

Wolfgang Propper has promised to establish a "Be- 
ginner's Column" but he has to apologize for this is- 
sue because of too many other engagements (lec- 
tures, courses). I think that Peter Mitic's comments on 
"Program-ming with DERIVE" might be useful even for 
not so experienced DERIVE users. Although Peter's 

article has "DERIVE 

3. x" in its title in my 
opinion it has not lost 
anything of its topical- 
ity even in times of 
DfW 4.09 and DfD 

4. x. 

I'd like to point out to 
you W. Weiskireh 's 
und H.Knechtei's re- 
quest for partners for 
an EU-COMENIUS- 
project. I am pleased 
to report that we have 
applied for a project 
"Modern 

Technologies Teaching Mathematics". With the Peda- 
gogical Institute of Lower Austria as the coordinating 
institution the following institutions have joined as pro- 
ject partners: ZKL Munster (Germany), CTM Univer- 
sity Plymouth (England), University of Gothenburg 
(Sweden), Freudenthal Institute (Netherlands), and 
the Catholic University of Leuven (Belgium). I'll keep 
you informed about the progress of the project 
(- if we will get the approval.) 

Finally I would like to call your attention to the en- 
closed programs and registration forms for the Get- 
tysburg Conference and for the TI-92 meeting in Mun- 
ster. 

With my best regards 


A bilef Histoiy of the DERIVE User Group 

From 1991 until December 1997 

1000 p^es - without the fmntp^es -containing: 

134 ooiiti'ibutio]i.s from 
^0 (lifrerent contributors from 
21 coimtries 

We had 

302 questions, answers, messages, comments^ 
j|]t,1he User Forum 

until DNLttZS we have put 

133 DERIVE- and Tl- related books on our bookshelf . 




Falls Sie Hire Mitgliedschaft noch nieht erneuert haben sollten, dann tun Sie das bitte so 
bald wie mdglich. Herzlichen Dank. 

If you have not renewed your membership for 1998 then please do that as soon 
as possible. Many thanks. 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. 

It must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the au- 
thor gives his consent for reprinting it in 
D-N-L. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Next issue: June 1998 

Deadline: 15 May 1998 


Preview: Contributions for the next issues 

Experimenting with GRAM & SCHMIDT, Schonefeld, USA 
Extracting Logic Propositions from Numerical Data, Etchells, UK 
3D-Geometry, Reichel, AUT 

Linear Programming, Various Projections, Word Problems, Bdhm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Fractals and other Graphics, Koth, AUT 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

Riemann, a package for the TI-92, Bohm & Propper, AUT/GER 

Parallel Curves, Wunderling, GER 

Quaternion Algebra, Sirota, RUS 

Concentric Curve Shading, Speck, NZL 

Drawing in the Plane with DERIVE, Mata & Torres, ESP 

Information Technologies in Geometry, Rakov & Gorokh, UKR 

Parametric 3D-Plots with DERIVE and 3DV, Welke, GER 

Various Training Programs for Students on the TI-92, Bdhm, AUT 

A Critical Comment on the "Delayed Assignation" := =, Kummel, GER 

and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Weth, GER; Kirmse, GER 
Aue, GER; Koller, AUT; Mitic, UK; Tortosa, ESP; Schorn, GER; 
Santonja, ESP; Dorner, USA and and and 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tent of 40 pages minimum. The goals of 
the DNL are to enable the exehange of ex- 
periences made with DERIVE as well as to 
create a group to diseuss the possibilities 
of new methodical and didactical ma nn ers 
in teaching mathematics. 

We have established a seetion dealing with 
the use of the TI-92 and we try to combine 
these modem teehnologies. 

Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)660 31 36 365 
email: nojo.boehm@pgv.at 
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HeinzRainer Gever, StKatharinen, Germany http://members.aol.com/gymgutenbg 

HeinzRainer.Geyer@t-online.de 

I produced another small tool, because I wanted my class to perform conversions using various num- 
ber bases. Unfortunately there are no functions HEX(a), OCT(a), DUAL(a) etc., which will convert a 
decimal number into another base. (I must admit, that I did not search very long). My results are no 
numbers, but lists of digits. Therefore it might be not very easy to continue with calculations. But it is 
just the beginning. Any suggestions are appreciated and welcome. 

#1: BASES. MTH 

#2: CaseMode Sensitive 

#3: REM(number, base) := MOD(number, base) 

SELCT(r) := [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, A, B, C, D, E, F] 

#4: r + 1 

#5: SEQU(number , base) := ITERATES(FLOOR(x , base), x, number) 

AUXCnumber, base) := VECTOR(SELCT(REM((SEQU(number , base)) , base)), k, 1, 

#6: k 

DIMENSION(SEQU(number, base)) - 2) 

CONVCa, x_) REVERSE_VECTOR(AUX(a , x_)) 

DUAL(a) := C0NV(a, 2) 

[DUAL(a) := C0NV(a, 2), QUAD(a) := C0NV(a, 4), PENT(a) := C0NV(a, 5)] 

[OCT(a) C0NV(a, 8), DEZ(a) C0NV(a, 10), HEX(a) C0NV(a, 16)] 

TESTS 

AUXC258, 16) = [2, 0, 1] 

HEXC258) = [1, 0, 2] 

HEXC253) = [F, D] 

0CTC27) = [3, 3] 

OCTC512) = [1, 0, 0, 0] 

DUALC27) = [1, 1, 0, 1, 1] 

C0NV(1199, 12) = [8, 3, B] 


As you can see there are some functions 
implemented in Tl-Nspire. Oh and Ob are the 
initials for hexadecimal and dual numbers. 



*Nicht gespeicherte 

258 ►Base 16 

Oh 102 


253 ►Base 16 

OhFD 


27^Base2 

ObllOll 


ObllOll^BaselO 

27 


[ 


4/99 1 


#7: 

#8: 

#9: 

#10 

#11 

#12 

#13 

#14 

#15 

#16 

#17 

#18 


There is no function for converting into any 
other number system. 


p 
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This is a DERIVE 6 program for converting decimal numbers into any other system. The re- 
sults are strings, of course. It should be a nice problem for students to program the conver- 
sion back to decimal numbers, Josef 

conv(x, b, r) := 

Prog 
r := "" 

Loop 

If X < b 

RETURN APPEND(IF(M0D(x, b) < 10, STRING(M0D(x , b)), 

C0DES_T0_NAME(55 + M0D(x, b))), r) 
r := APPEND(IF(M0D(x, b) < 10, STRING(M0D(x , b)), 

C0DES_T0_NAME(55 + M0D(x, b))), r) 

X := FLOORCx, b) 


[dual(x) conv(x, 2), oct(x) conv(x, 8), hex(x) conv(x, 16)] 

hex(258) = 102 

hex(253) = FD 

oct(512) = 1000 

dual (27) = non 

conv(1199, 12) = 83B 

conv(1199, 20) = 233 


Stefan Welke, Bonn, Germany Spwelke@aol.com 

A remark on John Alexiou's suggestion in DNL #25 concerning substitutions with the ITERATE func- 
tion. Alexiou mentions that his (very elegant) use of ITERATE improves the practice with the LIM- 
function insofar as it applies to any size and shape of vectors and matrices (especially) to number and 
symbolic calculations. I want to point out a special problem with symbolic substitutions which arises 
with simultaneous substitutions of several variables. 

I start with the definition of two functions: 

SUBSTl.MTH 

subst_a1ex(expr_, s_) := ITERATE(expr_, LHS(s_), RHS(s_), 1) 

subst_1 im(expr_, var_, s_) lim expr_ 

var_^s_ 

We now want to substitute in the expression y - x^ the variable x ^ a and y ^ b! " 

2 2 
subst_a1ex(y - x , [x = a, y = b]) = b - a 

2 2 
subst_1im(y - x , [x, y], [a, b]) = b - a 

We recognize that both functions are doing the same in this situation. The situation changes if the 
variable which should be substituted appears in the expression which substitutes the variable in ques- 
tion. Consider the following example: 
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2 2 
subst_alex(y -x,[x=x+y, y=x- y]) = - x + x - y 

2 2 2 
subst_lim(y - x , [x, y] , [x + y, x - y]) = - 4«x + x«(4«y + 1) - y - y 

We do not have the same result. The correct result should be: 

2 2 2 
(x - y) - (x + y) = - X + X.(l - 2-y) - y - y 

So neither SUBST ALEX nor SUBST LIM wins. Notice that both vesions work well if only one 
variable is substituted: 

2 2 
subst_a1ex(y -x , x=x+y)=-x 

2 2 2 
subst_lim(y -x,x, x+y)=-x - 2«x«y - y + y 

Conclusion: Simultaneous substitution works in both cases correctly only if the substituted variable 
does not appear in the expression which replaces this variable. 

DNL: Stefan Welke has submitted a contribution overcoming this difficuity to appear in the 
next DNL 

DERIVE 6 behaves as expected when applying its SUBST-command: 


2 2 2 
SUBSTCy - X , [x, y] , [x + y, x - y]) = - x + x-(l - 2-y) - y - y 

2 2 2 
SUBSTCy -x,x, x+y)=-x - 2*x*y - y + y 


The WITH-Operator of Tl-Nspire refuses performing the substitution. It needs auxiliary vari- 
ables: 


2 2 
y-x [x =a and y=b b^a ^ 


and y=x-y 

"Error: Invalid constraint" 

2 

a-b~‘\a=x-^y and b=x-y 

-x^+x{2y^l)^^+y 

I 

3/99 
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A.White. New York. USA 

Sirs, I am trying to calculate the double integral 


AWHITE@sbu.edu 


0.3^ 


0.37T 


j 5pe ^'dpds or equivalently ^5pe ^'dsAp. 


Q2k 


Doing these by hand yields approximatley 0.3038. Using DERIVE on a PC yields approximately 
0.3038. However, using Derive on a TI-92 yields 

INT (INT (5^p^#e^ (-p^s) ,s,2,inf) ,p, 0.2, 0.3) as a mess and 

INT (INT (5^p^#e^ (-p^s) ,p, 0 . 2, 0 . 3) , s, 2, inf ) as undef. 

Any suggestions on how to handle these on a Ti-92? 

Thank you. A1 

DNL: Compare your input with the input on the Tl-edit line. Both "INT" and "INF" are DERIVE 
- Syntax. Use the [oo] and the Integral key ([2nd] [7]). Take care of [e^], don't enter an ordi- 
nary "e". 


fFlT^Y fT? Y Y FHt Y fI Y Ffi ] 

f—|Rlgebra ICalc [other iPrgM 10 [Clear a-z... | 




I 1 in 


5 p e 

5 p e ^'^)cisdp 

5 p e '^'^)(isdp 


.3033 

undef 

.3033 





:s 1 

■j2j'2^5 p e ^ ^Jdpds 

.3G3771 



2 




■ ^ 


8*cosClln/120>#cos<29n/12O> 

MAIN RAD EXACT FUNC 1/ H 


The left screen is from a TI-92 (OS 2.07 of the V200 behaves the same). The right screen is 
from a Voage 200, OS 3.10. See below how Tl-Nspire copes with this problem: 


0.3 


( 

[5‘p e R^jdp ds 


0.2 


(0.740818)^-(,U+5.)(l. 10517)^-1. 5 (5+3.33333)j 


ds 


0.303771 


0.3 


ds 


0.2 


0.3 


0.2 


fs/^ e^^lds dp 


0.303771 


1/3 
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DERIVE 6: 


1 


0.3 - 2/5 - 3/5 

- p ■ s 5 ■ e 5 ■ e 

5-p-e dp ds = - 

0.2 2 2 


0.3 


1 


0.2 


- p-s 

5-p-e dp ds = 0.303771024S 


2 

0.3 

0.2 

0.3 

0.2 


- 3/5 1/5 

- p-s 5 - e - (e - 1) 

■ p-e ds dp = 

2 2 


J s.p.. 


1 


- p-s 

5-p-e ds dp = 0.303771024S 


Ivano Moschetti, Italy imoschet@novanet.it 

Hallo, Fm a beginner in DERIVE (and English). Em interested in algebraic expressions and other, that 
seems DERIVE bugs. What is the reason of this singular problem? 

For example: 

• X + 1 + — *— = 3 + —'—; DERIVE solve [x = 2] !?! Why?? 

x-2 x-2 

• ^ in graphic plot for x = 0, y isn't 1 !?! Why?? 

Inx 

Write to me other examples and reason. (My English is terrible, sorry) 

Ivano 

George F. Feissner, USA feissnerg@snycorva.cortland 

Moschetti's problem is an interesting one! When one solves the equation 
X + 1 + l/(x - 2) = 3 + l/(x - 2) 

using either the SOLVE command or the SOLVE function in DERIVE, one gets the spurious solution 

X = 2. 

I got the same result with 

x^2 + l/(x - 2) = 2 + X + l/(x - 2). 

Is DERIVE simplifying the equation by "cancelling" the l/(x-2) from each side? 
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A1 Rich, Honolulu, USA swh@aloha.com 

As X approaches 2, the left and right side of the above equations approach the same value. Since in 
DERIVE, 1/0 is well defined as ± oo, substituting x = 2 in the above equations yields equal infinite 
values on both sides of the equation. Thus, in DERIVE 's eyes at least, x = 2 is a valid solution. 

If you are unhappy with this situation, I recommend substituting solutions returned by DERIVE back 
into the original equation. If the result is infinite values on the left and right sides of the equation, you 
are free to disregard the solution as spurious. - Aloha, Albert D. Rich, Applied Logician 


DERIVE 6’s answers: 


SOLVE 


X + 1 + 


= 3 + — 

X - 2 X 




Cx = 2) 


SUBST 


LNCx) - 1 

j X , 0 I = 

L LN(x) J 


E 

3l 

■0.5 

-0.4 -0.2 

0.2 0 
■ -0 . 5 

-|— 1 Cross: 0 ^ 1 


Alfonso Poblacion, Valladolid, Spain alfonso@gauss.mat.eup.uva.es 

Problems with FOURIER 
Dear Derivers, 

Using DERIVE 3.0 I have found the following problem: 

I want to obtain The Fourier expansion of ABS(SIN(x)) in [-tt, tt] . Well, I loaded INT APPS.MTH 

2 5cost2xI 

and then simplified FOURIER (abs (sin (x) ) , x, -tt, tt, 2 ) . This gave me — - . If you 

K 3k 

plot this function together with the initial one, you can see at first sight that there is something that is 
strange in the behaviour of the second function. Then I simplified FOURIER (abs ( sin (x) ) ,x,- 

7 T , 71 , 4 ) and obtained — - + cos(3x) _ 2cos(4x) 

71 3n 4k 1571 

is not an expected FOURIER expansion. 


Then I tried to find the coefficients from its definition with INT (abs (sin (x) ) cos (nx) , x, - 

2 (- 1 )'^ 

TT, tt) /tt declaring n as a positive integer. Then you find — . If you build the expansion with 

7i(l - n) 

this result, this does not work either. (Not only seeing the plots, this gave me a first clue to suspect that 
something went wrong. By this time, I did the calculations myself by hand). 


Knowing the behaviour of the Fourier expansion, then I tried a third attempt. As the function is even, 
the coefficients can be obtained from 2/7r*int (abs (sin (x) ) cos (nx) , x, 0 , tt) . Then, finally 

2 ((- 1 )'^ + 1 ) 

you will have the right answer: — . 
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COS(0)/2 + COS(3*0)/6 + Cl = 0 , thus Cl = 1/3 

COS(0)/2 - COS(3*0)/6 + C2 = 0 , thus C2= -1/3 


Finally, we have: 




cos(x) cos(3x) 1 

2 6 "^3 

0 

cos(x) cos(3x) 1 


2 6 3 


for - 71 < X < 0 
for X = 0 
for 0 < X < 71 


The function I(x) is continuous and derivable on interval - tt <= x <= tt and, it is: dy/dx = I(x) 
So I(x) is the antiderivative function of y(x) on all interval [- tt, tt]. 

The Integral Fundamental Theorem gives: 


INT(y(x),x,-7r,7r) = I(7 t) -I(-7t) = -2/3 - (+2/3) = -4/3 that, of course, is correct. 

Now leFs see how DERIVE works. Unfortunately, DERIVE can not integrate functions defined as 
above. 

DERIVE uses, instead, same special functions as ABS(x), FLOOR(x), SIGN(x) etc. that makes the 
definition functions easier and more compact, but, on the other hand, hides same important conceptual 
problems. 

In fact DERIVE gives the following antiderivative function: 

INT (ABS (SIN (x) ) ^COS (2x) ,x) = SIGN (SIN (x) ) ^ (COS (x) /2-COS (3x) /6) 

We should believe that the function above is the antiderivative of y(x) on our interval. Well, we should 
get wrong! In fact the function SIGN(SIN(x))*(COS(x)/2-COS(3*x)/6) is not continuous (and , thus, 
not derivable) in x=0, tt ,- tt. We easy, can see it by the plot , as well, from the following: 

SIGN (SIN (tt) ) ^ (COS (tt) /2-COS (S^tt) /6) =” + -”1/3 
SIGN (SIN (0) ) ^ (COS (0) /2-COS (3^ (0) ) /6)=”+-”l/3 
SIGN (SIN (-7T) ) ^ (COS (-7T) /2-COS (3^ (-tt) ) /6) =” + -”1/3 

DERIVE uses limits to overcome the ambiguity of sign at the end of interval (and this is correct, 

I think) but can do nothing for the "internal jumps" as in 0. In fact if we name: 

DERIVE compute the following: 

LIM (SIGN (SIN (x) ) ^ (COS (x) /2-COS (3x) /6) , x, tt, -1) =-1/3 

LIM (SIGN (SIN (x) ) ^ (COS (x) /2-COS (3x) /6) , x, -tt, 1) =1/3 

INT (ABS (SIN (x) ) ^COS (2x) ,x, -TT, 7t)=-1/3 - ( + 1/3) = - 2/3 

I am not sure, but we should aspect wrong results when the antiderivative function makes "jumps" into 
integration interval. If it happens we must divide the original interval into sub-intervals (as we usually 
do by hand) and integrate separately. In our example, we have: 

INT (ABS (SIN (x) ) COS ( 2x) , x, -tt, Ti) = 

INT (ABS (SIN (x) ) COS (2x) ,x, -tt, 0) +INT (ABS (SIN (x) ) COS (2x) ,x,0,n) =-2/3- 

2/3=-4/3 

that, finally, is correct. 






D 


-N-L#29 


D E R I V E & TI-92 USER FORUM 


P 11 


There is another point: The function SIGN (SIN (x)) (COS (x) /2-COS (3x) /6) returned by 
DERIVE is an antiderivative of same function in -tt < x < ? 

The answer is yes if we introduce the generalized function delta(x) of Dirac: 

INT (delta (x) , X, -inf , inf ) = 1, delta(x) = 0 for x< 0 and x > 0 
With this statement, we can write: 

d/dx (SIGN (SIN (x) ) (COS (x) /2-COS (3x) /6) ) = ABS (SIN (x) ) COS (2x) + 2/3 
delta (x) ,for-7r<x<7i 

Integrating each member, for -7i< x < tt ,we have: 

d/dx (SIGN (SIN (x) ) (COS (x) /2-COS ( 3x) / 6 ) ) = 

LIM(SIGN (SIN (x) ) (COS (x) /2-COS (3x) /6) ,x, 71, -1) - 
LIM(SIGN (SIN (x) ) (COS (x) /2-COS (3x) /6) ,x, -7T, 1) = -2/3 

INT(2/3 delta (x) , x, -71, n) = INT(2/3 delta (x) , x, -co, co) = 2/3 
INT (ABS (SIN (x) ) COS (2x) , x, -71,71) = I 


So, rearranging we have: -2/3 =2/3 + 1 , thus I = -4/3 and that is correct 


The problem has been solved in the meanwhile as you can see below in the DERIVE 6 
expression: 

n 4 

/ |SINCx)| .C0SC2-X) dx = 

-Jl 3 


Finally, another aspect of built-in function SIGN(x) that makes confusion. For x=0 , it returns: 
SIGN(0)='’+-"1 . 1 hear some of my students saying: 

”the sign of zero is plus or minus one..!” 

And I have to scratch this silly idea by repeating continuously that "for each x there must be one and 
only one.etc., etc." 


Hope hearing other comment. Thanks . Hi, everyone. 
Leonardo 


Compare the DERIVE screen shots: DERIVE 4.x delivers the correct FOURIER expansion. 
Josef 



The respective DERIVE 6 screen shot is given on the next page. 
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(Sometimes you can follow an evolution of ideas in the Derive-News-Group. I find it useful for 
me - and I hope for you, too - to collect all the messages dealing with one problem and put 
them together. The following story is a good example of the fact that a question with four or 
five lines might bring a mathematical discussion into motion. Now you are invited to follow 


Wbt ^isitorp of Senrj) (Hranfjolm's; Ciiuation 


Tuesday, 20. January 1998, 16:28:44 henry.granholm@helsinki.fi 

The equation 4^ +7 = 2^“^ has one real solution (approximately - 0.868). 


If we solve the equation algebraically we get six complex solutions, but no real solution. There is 
nothing wrong in the complex solutions. We find all of them, and the real solution, if we use 
substitution y = 2^. (The equation is equivalent with (2"")^ + 7^2^ = 4, and the equation y^ + 7y = 4 has 
one solution in R^.) 

This problem was found in DfW 4.07 and 4.08. Can DfD find all solutions? Are there any other 
similar bugs in DERIVE? - Henry 

Tuesday, 20. January 1998, 18:40:19 jeff@friedmann.com 

I have run the equation on my HP 200 based DfD and gotten SEVEN solutions as indicated (not just 6 
complex ones), the first one given reduces to - 0.867943 on approximating (DERIVE solves the 
equation nicely without any numerical simplifying). - Jeffrey F.Friedman 


Wednesday, 21. January 1998, 18:05:41 Ja72@prism.gatech.edu 

From SOLVE(4^x+7 = 2^(2-x),x) with DfW 4.0 (and with DfD 4.06, too. Josef) 

we get: 


#9: 


X = -0.867943 
X = 1.43397 - 2.41271-i 
X = 1.43397 + 6.65200-i 
X = 1.43397 - 11.4774-i 
X = 1.43397 + 2.41271-i 
X = 1.43397 - 6.65200-i 


Approx ( #8 ) 


X 


1.43397 + 11.4774-i 


no BUG here - John 
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Wednesday, 21. January 1998, 17:49:51 henry.granholm@helsinki.fi 

feissnerg@snycorva.cortland.edu wrote 

> What exactly did you enter and what commands did you give? 

E.g. Solve - Algebraically on th eq. in question, Simplify-Basic on S0LVE(4^x + 7 = 2^(2 - x),x). 
lEs no difference. I get the exact solution on the complex roots, but nothing else. The result is the 
same if I try to solve (2'')^ + 7*2'' = 4. 

I donT get any error messages; nor any hints that there could be solutions DERIVE cannot find. 

After some testing (about a dozen equations) it seems as if the problem was connected with the 
equation y^+7y=4. The same problem seems to occur always when the equation has one positive irra- 
tional, and two complex solutions. There is no problem if the positive solution is rational, or, if the 
equation we get after the substitution has three real solutions. And there is no problem if the equation 
is of degree less than 3. 

Other problematic equations are e.g. -9" + 3^ ”"-2 = 0 

2-x_22x_4^0 

23x_2x^7 

Wednesday, 21. January 1998, 21:36:00 math@rhombus,be 

With DERIVE for DOS 3.1 the equation is solved symbolically into all 7 solutions, but the message 
"solutions not verified" appears! Apparently this bug creeped in from DfW 4.07 on. - Jan Vermeylen 

Wednesday, 21. January 1998, 23:07:33 Joe.H.Frisbee@usahq.unitedspacealliance.com 

A great deal more about Granholm's equation, 4" + 7 = 2^ ”", may be learned by separating the 
equation into its real and imaginary parts. In particular, letting the variable x be redefined as u + i v, 
both u and v real, permits equations for the two parts to be obtained and plotted separately. 

RE: 2^ (2u) ^COS (2^v^LN (2) ) -2^ (2-u) ^COS (v^LN (2) ) +7 = 0 

IM: 2^ (2u) ^SIN (2^v^LN (2) ) +2^ (2-u) ^SIN (v^LN (2) ) = 0 

The intersections of the two sets of solutions are the ones of interest for the original equation. It 
turns out that there are many more of these solutions than the numbers referred to in the previous 
messages. In fact, it appears that there is an infinite number of solutions. This is with the principal 
branch only and can easily be seen in the complex polar representation of y as well. 

It seems that the bigger question for any version of DERIVE is how it limits the number of 
solutions. By the way, I have used a DOS 3.03 and a WINDOWS 4.04 and have gotten from 3 to 9 
solutions depending on which form of x is used. All of my solutions are consistent with those 
previously reported. - Joe.H.Frisbee, Houston 
DNL: What a wonderful email address! 

Thursday, 22. January 1998, 09:17:15 swh@aloha.com 

When solving equations it is possible that spurious solutions may be introduced (for example, when 
both sides of an equation are squared). Therefore, after finding all candidate solutions, DERIVE back 
substitutes them back into the original equation, and only returns those that actually satisfy the 
equation. 
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Depending on the complexity of the solution and the equation, it can be very difficult to simplify the 
difference of the two sides of an equation to zero (this is called the zero-recognition problem in the 
computer algebra business). In your problem verifying the real solution to the equation requires the 
denesting of a cube-root of a sum of square-roots. 

A major redesign of the denesting algorithm implemented in version 4.07 of DERIVE significantly 
expanded the class of expressions that could be denested. However, a few expressions formerly 
denestable were no longer. 

Unfortunately, the verification of the real solution of your equation results in one of these expressions. 

The problem will be resolved by version 4.09 of DERIVE. A DERIVE for Windows version 4.09 
update will be available for downloading from www.derive.com as soon as it is released. 

We appreciate your bringing this problem to our attention so it could be resolved. 

>Can Derive for DOS find all solutions? 

Since any given version of DERIVE for Windows and DERIVE for DOS use the same math engine, 
they will both return the same results for a math problem. 

>Are there any other similar bugs in Derive? 

The zero-recognition has been proven to be insoluable in general. Thus, there will always be problems 
that any computer algebra system, including DERIVE, cannot do. - Aloha, Albert D. Rich 

(DERIVE 4.09 shows seven solutions, indeed. Josef) 

Thursday, 22. January 1998, 10:57:20 J.Wiesenbauer@tuwien.ac.at 

To All, as the flood of e-mails keeps pouring in concerning Henry Granholm's equation 

4^x+7 = 2^ (2-x) (^) 

I strongly feel that I should also make a contribution to it. As Henry has already pointed out one could 
solve that equation by substituting y = 2^x and solving the resulting equation 

y^2+7 = 4/y 

Subsequently for each of the three different solutions y = yl,y2,y3 of the latter equation the equation 

2^x = y 

must be solved. If yl is the unique real solution of (**), then 

xl=LN( (SQRT (1353) /9 + 2)^(l/3) - ( SQRT ( 1353 ) / 9 - 2) ^ (1/3) ) /LN (2) 
is the corresponding unique real solution of (*). Surprisingly enough, after substituting x = xl in (*) 
you will get the result 'wrong', i.e. according to A1 Rich DERIVE fails to solve the corresponding 
zero-recognition problem. 

This accounts for the disappearance of the real solution. But what about the complex solutions? After 
all, the corresponding zero-recognition problem is a lot more complicated for them and again you can 
show by resubstituting the solutions x = x2 and x = x3 that DERIVE fails to solve it. This is what I 
think might be a plausible explanation: In contrast to the real solution DERIVE doesn't even try to 
verify them! 

Strangely enough, there is a far more serious bug in connection with this problem that has been 
overlooked by all contributors (with the possible exception of Joe Frisbee). Which one? Well, 
obviously the exact solution of (***) w.r.t. the variable x should be 
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X = log (y, 2 ) +2k'^pi'^#i/ln (2 ) , 

where k is an arbitrary integer. For some reason DERIVE ignores all solutions except for the one with 
k = 0, if log(y,2) is real, and k = -1,0,1 otherwise. Why? 

I haven't the faintest idea! Only God and A1 Rich know the answer... At any rate, this accounts for the 
6 complex solutions in the case above. (It goes without saying that I am using the current version 4.08 
ofDfW.) 

Let me conclude by pointing out that I enjoy problems of this kind a lot (thanks to Henry!) and I 
would like to see more of them here! - Cheers, Johann 


Friday, 23. January 1998 00:06:02 Joe.H.Frisbee@usahq.unitedspacealliance.com 

Following Johann Wiesenbauer's lead, I have formalized what I believe to be a complete and 
exact solution description. There are two families of solutions, one in the left half plane (LHP) and the 
other in the right half plane (RHP). The previously presented real solution (- -0.8679) is a member of 
the LHP family and the reported complex solutions are in the RHP family. The LHP family is 
expressed rather simply and it is odd that DERIVE does not give additional solutions in this family. 
Using the complex representation, X = U + V i , the two families of solutions (provided I have not 
mistyped) are: 

U_LHP = (2*LN((SQRT(451) + 6*SQRT(3))^(l/3) - (SQRT(451) - 6*SQRT(3))^(l/3)) - 
- LN(3))/(2*LN(2)) 

VLHP = (2*N*PI)/LN(2) 

U RHP = (LN((559 - 12*SQRT(1353))^(l/3) + (559 + 12*SQRT(1353))^(l/3) + 7)- 
-LN(3))/(2*LN(2)) 

V RHP = PI*(2*N + 1)/(2*LN(2)) + (ASIN(1 - 3*U_RHP)*(-1)^N)/LN(2) 

Notes: 1) In these equations N is an integer on (-INF,+INF). 

2) The real component of X is constant in each family. 

3) The arcsine function presents the solution on (0,PI/2). 

4) N = 0 yields the single real solution which is in the LHP family. 

5) U LHP is the same as Johann's XI solution of (*). 

I would like to echo Johann's comment concerning interest and challenge in such problems. As a 
user of DERIVE, I seldom concern myself with its inner workings. Problems like this should help to 
keep us all aware of the complexities of mathematical problems and the software capabilities required 
to handle them. Also, I would like to say that we have not really been talking about a bug in DERIVE 
(apologies to Soft Warehouse) but more of a level of sophistication. I am continually pleased and 
amazed at what DERIVE can do. It is an exceptionally useful tool for applied mathematical analysis, 
as well as for education in which many of you are involved. 
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Friday, 23. January 1998 17:50:30 J.Wiesenbauer@tuwien.ac.at 

To All, 

I apologize to all who are sick and tired of getting e-mails concerning Henry Granholm's notorious 
equation 

X 2 - X 

4 +7 = 2 

but after reading Joe Frisbee's last e-mail I have a terribly guilty conscience and I simply want to get 
rid of it. 

It appears that there is something wrong with Joe's solution of the equation above. The thing is that his 
second series U RHP + #i*V_RHP doesn't make sense since V RHP, which is supposed to be real, 
yields complex values due to the fact that the argument of the arcsine function is out of the range [- 
1,1]. Since he is referring to my "lead", this leads to the most embarassing conclusion that there must 
be something wrong with it. 

To make up for it, I have written a DERIVE-routine for equations that can be transformed by a 
substitution of the form y=a^x (assuming that x is the main variable) into a solvable one: 


EXPSOLVE (u,x,a) 
Check e.g. 


VECTOR (x = LOG(y_,a) + 2-PI-0-#i/LN (a) , y_, 
RHS (SOLVE (LIM(u,LOG(y_,a) ) , y_) ) ) 


EXPSOLVE (4^x + 7 = 2^(2 - x), x, 2) 


as well as Henry's "other problematic equations": 


EXPSOLVE (- 9^x + 3^(5 - x) -2=0, x, 3) 


EXPSOLVE (2^ (3x) - 2^x = 7, x, 2) 

EXPSOLVE (2^ (-X) - 2^(2x) - 4 = 0, x, 2) 

(Note: If you want to see the approximate solutions after simplifying the expressions above you should 
press the "Approximate"-icon twice. This will lead to a better representation!) 

But the real value of the routine above should be that you can use it as a "tutorial", as it were, on 
solving such equations manually by studying the code! 


For the rest, I agree wholeheartedly with Joe's opinion that the term "bug" might be inappropriate in 
the case in question. Frankly, I didn't expect DERIVE (or any other CAS for that matter) to solve such 
equations correctly... 

Ah, I feel much better now! Apologies again... - Cheers, Johann 


Friday, 23. January 1998 19:53:46 Joe.H.Frisbee@usahq.unitedspacealliance 

To All, 

Well I warned you there could be typing errors. Johann you are correct that the V RHP is in error. 
The arcsine term should read ASIN(2^(1-3*U_RHP)) which will satisfy the [-1,1] requirement. With 
this correction the right half plane V is: 

V RHP = PP(2*N + 1)/(2*LN(2)) + (ASIN(2^(1 - 3*U_RHP))*(-1)^N)/LN(2) 

Hopefully this will put everything in a complete and errorless form. Now if I could only find some 
software that compares my intent to my action I'd really be happy! 

Regards from: Joe Frisbee 
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In the following you can see how DERIVE 6.10 performs now solving Granholm’s Equation. 
The solutions in Approximate Mode are very strage! 


Precision := BB BBS8HBEEB 

Notation := Decimal 


-0.S6 79410104 


X 2 - 

S0LUTI0NSC4 +7=2 



CO + 9.0647202S3.L 
CO - 9.0647202S3.L 
CO + 13.5970S042.i_ 

CO + 4. 532360141- i. 

CO - 4. 532360141- i. _ 


Now I am working in Exact Mode. The SOLUTIONS matrix is the result of approximating the 
exact results - which are looking very “exotic” 


Exact 


Notation := Rational 

r X ^ ^ 1 

L50LUTI0N5C4 +7=2 , x)J^ 

-0.S679410104 

1.433970505 + 2 . 41271701S- i. 

1.433970505 - 6 . 652003264 ■ i. 

1.433970505 + 11 . 47743730 ■ i. 

1.433970505 - 2 . 41271701S- i. 

1.433970505 + 6 . 652003264 ■ i. 

1.433970505 - 11 . 47743730 - i. 
-0.S679410104 + 9 . 0647202S3 ■ l 
_ -0.S679410104 - 9 . 0647202S3 ■ l _ 


Josef Lechner, Viehdorf, Austria 


^ I — 

r '] 

Algebra 

r 1 

Calc 

r FH’r 1 

other 

r i 

PrgnlO 

f Ffi 1 

Clear a-z... 


■ taylorlf ^ , x , I gJ 

x^*^ x^ X® > 

3715S912Q0 1S579456G 1G32192G 64^ 

T x^*^ x^ X® 

■J 1. 37158912GG 18579456G 1G32192G i 

^ -2 ► 


If <ans<l>,x,B,l> 


MAIN RflD EKflCT FUNC E/SO 


iosef.lechner@online.edvg.co.at 

A bug on the TI-92??? At least a strange 
behaviour! ! 

the indefinite integral: no problem, but: 
the definite integral: busy, busy, busy, 


See more on page 43. 
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Betrachtung des relativen Fehlers bei einer Standardaufgabe zum 
zweiten Strahlensatz mit Hilfe des TI-92 

Inspection of the relative Error in a Standard Problem to the 
2nd Theorem of Rays using the TI-92 

Jorg Mali, Bonn, Germany 


Eine bekannte Aufgabe zum 2. Strahlensatz lautet: 


"Eine Erbse von 6mm Durchmesser verdeckt 
gerade den Vollmond, wenn man sie 66cm vom 
Auge entfernt halt. Berechne einen Schatzwert 
fur den Mondradius." (LS 9, Klett Verlag, S 145). 

"A pea with 6mm diameter hides the full moon if 
you hold it at a distance of 66cm from your eye. 
Find an estimate for the moon 's radius. " 

Aus dem zweiten Strahlensatz erhalt man einen Schatzwert von ca 1745.5 km. Vergleicht 
man diesen Wert mit einer Lexikon-Angabe fur den Mondradius von 1738 km (Meyers 
Grofles Taschenlexikon, Mannheim 1987), so errechnet man einen relativen Fehler von 



1745.5-1738 

1738 


100 % 


0.43%. 


The 2nd Theorem of Rays leads to an estimate of 1745.5 km. In a lexicon we find 1738 km 
for the moon's radius, so we have a relative error of 0.43%. The statement "66 cm" could 
puzzle you. There is somebody checking if the pea is reaiiy hiding the moon and 
simuitaneously trying to measure the pea 's distance from his eye. So the question appears 
how the estimates' errors would change if somebody wouid (mistakeniy) beiieve that the pea 
is hiding the moon at 65 cm or at 68 cm. 1/1/e define together with the students a "Moon 
Estimation Function" y1(x) with x the distance eye - pea. The table on the TI-92 shows that 
this function is monotonous - as expected - and that the error depends on over- or 
underestimation of the pea 's distance. 


Die Angabe "66 cm" konnte stutzig machen: Da versucht jemand zu erkennen, ob eine Erbse 
nun wirklich gerade den Mond verdeckt und unternimmt gleichzeitig MeBversuche mit einem 
Zollstock vor seinem Auge. Es liegt also die Frage nahe, wie sich die Fehler der Schatzwerte 
verandern, wenn man (irrtumlich?) glaubt, dass die Erbse den Mond schon bei 65 cm oder 
bei erst 68 cm verdeckt. Dazu kann man mit den Schulern eine "Mondschatzfunktion" 
definieren, die einem Augenabstand x [in cm] einen Schatzwert fur den Mondradius y1(x) [in 
km] zuordnet: 


384000 


• 0.000003 


115200 


100000 
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Ein Blick auf die Wertetabelle des TI-Q2 zeigt 
das erwartete Monotonieverhalten und dass 
es fur den absoluten Fehler bedeutsam ist, 
ob man den Augenabstand uber- Oder 
unterschatzt. 


■FIT^HOT Y : S T-"" 

T Zoon 

Y Y 

: s- 


*PLDTS 

. , 30400G , ^ 

■^yl= 3. e-6 

X 

yl 

64.000 

1300.0 

65.000 

1772.3 

100000 

y2= 

y4= 

y5= 

y|= 

y7= 

66.000 

1745.5 

67.000 

1719.4 

63.000 

1694. 1 

69.000 

1669.6 

70.000 

1645.7 



wl GG/<X/1GGGGG>*G - GGGGG3 

INTEG RAD AUTD FUNC 


Die "Fehlerprozentfunktion" kann man mit den Schulern unter Ru ckg riff auf y1 entwickein: 


y2(x) = 


|y1 (x)-1738 
1738 


• 100 


Der Graph dieser Funktion ist interessant: 


We produce an "Error Procent Function" using y1 and study its interesting graph. We can 
find some connections between mathematical facts and reality, eg 

- Can you explain the asymmetric behaviour? Does it follow your experience? 

- How can you interpret a point (xc\yc) presented in the "Trace"-mode? 

- What is the solution of y1(x) <5? How can you find the solution graphically using the 
TI-92? 


Auch hier sind Verknupfungen von mathematischen Sachverhalten mit der Wirklichkeit ge- 
geben, beispielsweise: 



xmin= - 5 . 

xmax=lO0. 

xscl=10. 

ymin= "5. 

ymax=100. 

yscl=10. 

xres=5. 


\ 

^ 



MAIN RA 


KACT FUNC 


-Was bedeutet die Asymetrie? Foigt sie aus der An- 
schauung? 

-Wie kann man einen Punkt (xc|yc), der im "Trace" 
Modus angezeigt wird, interpretieren? 

-Was gibt die Losungsmenge von y2(x) < 5 an? Wie 
kann man sie mit Hilfe desT/-92 ermittein? 


Schliedlich kann man den TI-92 das Minimum der Funktion numerisch ausrechnen lessen. 
Es ergibt sich, unabhangig vom Untersuchungsintervall etwa T(66,28308|0, 00000653). Wie 
ist dieser Punkt inhaltlich zu deuten? Ist er wirklich ein Tiefpunkt? Eigentlich muUte doch - 
anschaulich gedacht - das Minimum 0 sein. Das sieht man auch auf einer anderen 
Darstellungsebene, namlich anhand der Formel: 

|y1(x)-1738| 57600 

J ^-100 = 0 « y1(x) = 1738 ("Solve"-Befehl) x = 

1738 869 

Die Schuler erkennen hier nicht nur die Fehibarkeit 
des Rechners, sondern machen die weitgehende 
Erfahrung, dass man doch ein exaktes Ergebnis 
erhalten kann, wenn man dem Computer die rich- 
tigen Fragen stellt. 

(Finding the zero(s) of y2(x) without the abs- 
function could be another approach. Josef) 


fFiT^Y Y Y Y Y Y*"^ .^Y 1 

T Zoom Trace ReGraph Math Draw ^ ^ 

■ soloe(y2(x) = 0 , x) 

_ 57600 

y ' ' ' Y 

Zero -X. 

xc: 66. 233034 ... 

369 

■ soloe(y2(x) = G , x) 

X = 66.2330340 


yc: 0. 0000000 . 

INTEG RAD EKACT FUNC I 









Finally the TI-92 gives the numerical value of the function's minimum 
(66. 28308\0. 00000653). Can you explain that point? Is that really a mimimum? It should have 
the value OH Using an equation we find our consideration confirmed. We obtain an exact so- 
lution. So the students do not only recognize the fallibility of the machine but also the fact 
that it might be necessary to put the question in the appropriate form. 


You can do that all using DERIVE, too. 

Produce your own TABLE-function (tblset in- 
cluded); 

I'd like to add another question to the stu- 
dents: 

Can you explain the endbehaviour of 
y2(x)? 

COMMttMD: Center Delete Help Meue Options Plot Quit Range Scale Transfer 

Uindou aXes Zoom 
Tracing expression #6 

Cross x: 66. 4062 y:0 Scale x:50 y:50 Deriue 2D-plot 

The DERIVE Table function: 

TABLE (u_, x_, xl_, x2_, dx) : =VECT0R ( [x_, u_] , x_, xl_, x2_, dx) 

EXPLOITING NEW FEATURES IN DERIVE 3: MULTIPLE DECISIONS 
AND WHOLE STRUCTURE PROGRAMMING 

Peter Mitic, Southampton, UK 

ABSTRACT 

Two important features were introduced in DERIVE, version 3, which simplify and extend 
programming considerable. The first is the provision for a multiple branch decision construct. 
The second is a new form of the vector construct. The latter allows a whole structure 
programming paradigm, and both can simplify programming considerably and produce faster 
code. An example from data analysis is used to illustrate the effectiveness of these 
constructs. 


INTRODUCTION 

Programming in DERIVE is often limited by a lack of programming structures which are 
present in many other procedural programming languages and symbolic computation 
packages. We suggested the addition of several features to the DERIVE language, and, 
following some discussion with Soft Warehouse (Al Rich 1994), two specific suggestions 
have now been adopted in version 3. Specifically, we requested a multiple decision 
construct, similar to CASE in Pascal, and a loop construct such as FOR in C. The 
equivalents of CASE and FOR are now possible, albeit indirectly. 
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We show how these constructs can simplify some programming tasks considerably. Two 
examples are presented to illustrate them. The first is to generate a list of random numbers 
which does not have a uniform distribution. This process has a clear educational advantage 
in that it teaches a specific technique for transforming uniformly distributed random numbers 
in a simple way. The second is to produce a histogram of the data produced. Both are set in 
the context of simulating random events, such as throwing a die. DERIVE is also compared 
with Minitab in this respect. 


SIMULATION OF RANDOM EXPERIMENTS 

A random experiment is a trial for which there is a well defined set of outcomes, but the result 
of any particular trial cannot be determined in advance. Examples are throwing a die, tossing 
a coin, noting the manufacturer of the next vehicle which appears at a particular point in the 
road, or measuring the height of an object. A random variable is associated with each 
experiment. The internal DERIVE function we use for simulating random events is random, 
which may be adapted to simulate discrete or continuous events, random(z) + l generates a 
random integer between 1 and z inclusive and random(I) generates a random seal number 
between 0 and 1 . Both are uniform distributions. The following session shows how a 
sequence of 10 dice throws and a sequence of 10 random numbers in the range (0,1) may 
be simulated. 

#1: RANCOHCO) = 1543547412 

#2: die6 ;= VECTOR C RAN DOM (6) + 1, i, 1, 10) 

#3: die6 = [5, 2, 4, 2, 1, 1, 5, 2, 5, 3] 

This may be adapted to generate random 
numbers in the range (a,b) by replacing 
RANDOM(l) bya + (b-a)*RANDOM(l). 

Random number generation by DERIVE is 
undoubtely slower than the random number 
generation by a compiled or dedicated 
package, but is convenient to use this facility 
in conjunction with other algebraic facilities, 
and to avoid tedious import and export of 
data. 


#4: 

#5: 


# 6 : 


urn fo rmOl ;= VECTOR (RANDOM (1), 1, 10) 

[unifbrmOl] ’ 

0.1720B374S0 
0.302S220951 
0.3251B5360S 
0.7732B77516 
0.1S9SB34924 
0.41210S0S49 
0.733257S500 
0.35357024S7 
0.6601174S69 
0. 999125 S9SS 


Don't forget to start each session with simplifying RANDOM (0) to make your random numbers really 
"random". Otherwise you would get the same sequence of random numbers each session. 
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THE WHOLE STRUCTURE AND RECONSTRUCTION PROGRAMMING PARADIGMS 
FOR TRANSFORMING DATA 

Given a list, suppose that we wish to transform every element in the list in the same way. For 
example, we may wish to square each element. The pseudocode for this process is as 
follows. 


Procedure Transf orm ( v) 

For i:= 1 to n 

Transform element i of vector v 
EndFor 
EndProcedure 

Transforming every element in a vector was a tedious process in DERIVE, version 2, and 
involved isolating each element in the vector by deconstruction, transforming it, and then 
reconstructing the vector. A template for this reconstruction programming paradigm is: 

VECTOR(F(ELEMENT(v,i_)) ,i-,l,DIMENSION(v)) , 

where v is the vector to be transformed, and F() is a transformation function such as 
F(x):=x'^2. 

The improvement we suggested is to introduce a new programming paradigm which avoids 
deconstruction, transformation and reconstruction. In this paradigm, the vector to be 
transformed and the transformation function are merely combined in a single function. The 
syntax which has been implemented is vector(f(x) ,x,v) . This has two principal 
advantages. First, it greatly simplifies programming syntax. Secondly, it works at least three 
times as fast as reconstruction. 

Given that the numbers generated by DERIVES RANDOM function are uniformly distributed, 
we can transform a list of such numbers such that the resulting numbers have a non-uniform 
distribution. In the case where the distribution function, F(x), for a random variable, Y, which 
has a non-uniform distribution, is one to one, it is easy to transform observations from a 
uniform distribution described by a random variable X so that they have the distribution 
function F(x). A proof of this result may be found in, for example, (Tocher 1963) and (Mitrani 

ll?§^ example below we generate random numbers which are exponentially distributed with 

In (1 - X) 

probability density function f(X) = 2 e'^^ (X> 0). The required transformation is Y = — . 

The reconstruction paradigm is illustrated in expressions #10 - #12, below, and the whole 
structure programming paradigm is illustrated in expressions #14 - #17. 

Loca - x) 

# 7 : FCx) ;= 

-2 

#S: urn fo rm03 ;= VECTOR CRANCOH a), i, 1, 5) 

#9: umform03 ;= [0,0376575, 0.S3266, 0.3571A5, 0.S07229, 0.A5A272] 

expranB := VECTOR CF(unifbrm03 ), i, 1, DIM(unifbrm03)) 

# 10 : 1 

#11: expranB 

#12: [0.0191924, 0.S93S63, 0.22091S, 0.S23126, 0.302S17] 
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#13: 

aa 




#14: 

expran33 := VECTOR(FCx), x, unifbrm03) 


#15: 

expran33 




#16: 

[0.01919247973, 

0 . S93S640492 , 

0.22091S4065, 0 . S2312S33S9 , 

0.3028175 5 SI] 

#17: 

[0.01919247973, 

0 . S93S640492 , 

0.22091S4065, 0 . S2312S33S9 , 

0.3028175581] 

#1S: 

uniform05 := VECTOR (RANDOM (1) . 

1. 1. 5) 


#19: 

expran := VECTOR(F(unifbrm05 ) 

■i 

,1,1, DIH(unifbrm05)) 



#20: expran 

#21: [0.30S353, 0.26S722, 0.292794, 0.768173, 0.461253] 

#22: expranl VECTORCFCx), x, utrifiormOS) 

#23: exp rani 

#24: [0.4225S4, 0.604105, 0.509711, 0.277994, 1.12106] 

#25: [0.662413, 0.69277, 0.2S5336, 0.10150S, 0.3S7406] 

#26: NotationDigits := 6 

Do you note the difference: we have the same random numbers from #12 to #17 (because uni form03 
was first evaluated and the generated numbers then used several times, uni formOS (#18) was not 
evaluated and then each call of uni form05 in other expressions (here expranl) result in a new 
sequence of random numbers. Josef 


The compactness of the new programming construct is readily apparent. The new vector 
construct works between three and six times as fast as the deconstruction-transformation- 
reconstruction technique because considering each element in a list and transforming it is 
■hard wired' in the DERIVE.EXE file. 

More importantly, this is the first application in DERIVE of the concept that it is better to treat 
an object as a whole rather than consider individual elements in the list. This concept has 
been recognised for some time in Maple (with the function map) ans is a major design 
feature of Mathematica (with the function Map). This is, however, only a partial solution to the 
'loop' problem, since it is now possible to use the new vector construct to implement a loop 
in DERIVE if and only if all the elements within the scope of the loop are contained in a single 
list. 

Any suitable package can, of course, be used to generate random numbers from a given 
(non-uniform) distribution. Minitab, for example, has a wide range of options for common 
distributions, and is only necessary to choose one and generate as much data as is required. 
Minitab is also usually faster than DERIVE. However, this does not emphasise the theoretical 
considerations behind generating the random number distributions concerned, and we 
wished to stress the relevant concepts for later work. 
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THE MULTIPLE DECISION CONSTRUCT AND FREQUENCY COUNTS 

It proved to be very difficult to produce grouped frequency counts for a vector of random 
data, as generated above, despite the simplicity of the problem. In a simple example 
involving discrete data, given a sequence of throws of a 6-sided die, we wish to count the 
number of ones, twos, threes etc. For continuous data, we might wish to assign random 
numbers in the range (0,1) to classes (termed buckets) 0 - 0.25, 0.25 - 0.5, 0.5 - 0.75 and 
0.75 - 1. This was such a problem in DERIVE, version 2, that for 100 random numbers or 
fewer it was preferable to produce grouped frequency counts by hand. In order to program 
this, we would wish to implement pseudo-code such as following, in which dice throws are 
counted. 

Procedure CountVector (v) 

For i:= 1 to 6 

Initialise Count (i) to zero 
EndFor 

For i:= 1 to Elements_in_vector (v) 

Case Elements_in_vector (v) is 
1: Increment Count (1) by 1 
2: Increment Count (2) by 1 
3: Increment Count (3) by 1 
4: Increment Count (4) by 1 
5: Increment Count (5) by 1 
6: Increment Count (6) by 1 
EndCase 
EndFor 
EndProcedure 

The problem is that it was not possible to implement a construct similar to the CASE 
statement in Pascal directly. In addition, DERIVE does not support external memory 
variables such as Count in the pseudo-cose above. Consequently, such memory variables 
must be maintained internally to a function. A multiple IF statement can be written, but it 
rapidly becomes very unwieldy and difficult to alter. This makes programming difficult and the 
result is slow to run. One implementation for version 2 is given in Appendix 1, which is a 
listing of a suitable utility file, COUNT.MTH. This contains a function POS(e,v) , which 
locates the position of an element e in a vector v. This function is an important function in list 
processing because it can be used extensively for pattern matching. 

In the session below, we simulate 10000 throws of a die. The bucket boundaries are stated 
as [1 .5, 2.5, 3.5, 4.5, 5.5] and the vector of random numbers can be counted directly. We 
also generate 10000 uniformly distributed random numbers in the range (0,1) and classify 
them into buckets 0-0.25, 0.25 - 0.5, 0.5 - 0.75 and 0.75 - 1. The vector buckets is 
redefined. The boundaries 0 1 and 1 are implicit in that a number which is less than 0.25 
necessarily correspons to the bucket 0-0.25 and a number which is greater than 0.75 
necessarily corresponds to the bucket 0.75-1. The function POSITIONS must be called in 
conjunction with count. 


You have to preload the Utility file COUNT.MTH (Appendix 1). 
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(The sample with 10000 throws was simulated by me using DfD 3.1 4XM and DfD 4. The 
calculation times shown in the comments are from versions 3.14 and 4, and DERIVE 6, 
resepctively. Josef) 


#1: LOADCD: \D0KUS\DNLs\DNL9S\HTH29\C0UNT.HTH) 

#2: bucketsl := [1.5, 2,5, 3.5, 4,5, 5,5] 

#3: COUNT CVECT0RCRAND0MC6) +1, -i , 1, 10000), bucketsl) 

#4: [1678, 1630, 1691, 1671, 1683, 1647] 

27.5 sec DfD 3.14, 34.8 sec DfD 4, 1 .1 7 sec DERIVE 6 
#5: buckets2 := [0.25, 0.5, 0,75] 

#6: COUNT (POSITIONS CVECTOR (RANDOM a), i, 10000), buckets2) , buckets2) 

#7: [2512, 2515, 2511, 2462] 

In DERIVE, version 3, it is now possible to write a multiple IF statement in a compact form, 
although the CountVector pseudo-code still cannot be implemented because memory 
variables outside procedures are still not supported. The functions positions and count, 
listed in Appendix 1 , can be combined in the following definition. 

POSiTiONSlO then simplifies to the same type of vector as expression #10, above. 
POSiTiONSi is a less flexible expression than the combination of positions and count, 
because it is hard to alter the bucket boundaries by merely stating new ones. However, this 
amendment permits generation and counting of random numbers in about one fifth the time 
that was taken previously. It is therefore worthwhile. 

#8: POSITIONSl(v) I([IF(x < 1.5), IF(1.5 < x < 2.5), IF(2.5 < x < 3.5), 

IF(3.5 < X < 4.5), IF(4.5 < x < 5.5), IF(5.5 < x)], x, v) 


#9: P0SITI0NS1(VECT0R(RAND0M(6) + 1, i, 1, 10000)) 

#10: [1682, 1708, 1707, 1628, 1629, 1646] 

12.8 sec DfD 3.14XM, 13.4 sec DfD 4, 1.75 sec DERIVE 6 

In order to increase the usefulness of positionsi, it is necessary to note that construct 
such as IF (alpha < x < beta) occur repeatedly in POSITIONSI. This string is 
prepended and appended by terms iF(x<alpha) and if (beta<=x) respectively. In the 
following function, the bucket boundaries are held in a vector b and vector is used to 
implement the repeated string of IF constructs. 

P0SITI0NS2(v, b) ;= ICriF(x < b ), VECT0R(IF(b < x < b ), f, 1, DIM(b) - 1), 

1 1 1+1 


# 11 : 
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The result is a nested vector of the form [ei, [02, 03, ... en-i], On], which is iess usefui than the 
unnested form [ei, 02, 03, ... en-i, On]. This nested form may be 'flattened' by applying the 
function flatten, which removes the inner brackets by appending the third and prepending 
the first element of v to its second element. 

#12: P0SITI0NS2(VECT0R(RAND0M(6) +1, f, 1, 10000), bucketsl) 

#13: [1731, [1686, 1652, 1620, 1666], 1645] 

FLATTEN(v) := APPEND([v 1, APPEND(v , [v 1)) 

#14: L iJ 2 L 3j 

#15: FLATTEN ([1731, [1686, 1652, 1620, 1666], 1645]) 

#16: [1731, 1686, 1652, 1620, 1666, 1645] 

Combining positions2 and flatten, we can define positions as below: 

#17: P0SITI0NS3(v, b) := FLATTEN(P0SITI0NS2(v , b)) 

#18: POSITIONS3(VECTOR(RANDOM(6) +1, 1, 1, 10000), bucketsl) 

#19: [1636, 1625, 1707, 1709, 1642, 1681] 

22.8 sec DfD 3.14XM, 35.1 sec DfD 4, 2.95 sec DERIVE 6 

The price for using the new function positions3 is that it takes about twice as long to 
execute as the less flexible positionsi. 


ILLUSTRATIONS FOR FREQUENCY DISTRIBUTIONS 

In order to proceed with further work on data analysis, we required a means to graph a 
histogram of data obtained as above. The code for producing a rudimentary histogram is 
discussed here for the completeness. 

Many packages, including the majority of dedicated statistics packages and spreadsheets, 
and some symbolic manipulators, have built-in functions for graphing and producing 
histograms. DERIVE does not, and specific functions have to be coded. The result can be 
slow and inflexible in that non-equal interval widths are harder to code. However, the equal 
interval widths option is standard in many other packages, and anything more sophisticated 
is not. 

Functions for producing a histogram in which the interval widths are constant are contained 
in BARCHART.MTH, listed in Appendix 2 . We illustrate the case of graphing the grouped 
frequency distribution of 10000 random numbers which simulate tossing a fair die 10000 
times. The output (# 13 ) obtained by simplifying 

VBARS ([1655 ], buckets, 0 . 8 ) is not given because it is long and only serves to provide 

data for producing the graph. 


#20 : LOADCD : \D0KUS\DNLs\DNL98\MTH29\BARCHART . MTH) 

#21: VBARS([1636, 1625, 1707, 1709, 1642, 1681], bucketsl, 0.8) 
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DERIVE for DOS 


1400 



DERIVE 6.10 


Alternatively, the problem can be solved by exporting a vector as a text file, manipulating it 
using an external program and then importing the result. We have used the same method 
(Mitic and Thomas 1995) to remove unwanted symbols from a vector. This was succesful, 
but at the expense of a non-negligible Pascal programming overhead. In this case, it would 
only be worth pursuing if the time gained from using a compiled language is greater than the 
time lost by creating the export text file and using it elsewhere. DERIVE does have the 
advantage that it provides an integrated environment in which simulation, graphing and 
algebraic manipulation can all be done. 


CONCLUSION 

Two new advances to the DERIVE programming language have been identified and 
examples of their efficacy have been presented. Although they are not a complete solution to 
the loop construct and to the multiple decision construct which we originally envisaged, they 
can be applied with reasonable ease, and adapted to more general cases. We intend to use 
the simulation and illustration activities presented here to illustrate some concepts in 
probability and statistical models in a later paper. 

Using DERIVE is not without its problems and in some cases, computations are better done 
by hand, usually because this is quicker. This is clearly not practical with a large amount of 
data, so the need to solve problems by programming becomes important. Random number 
generation, classification of data and graphing can all be performed faster and with simpler 
programming by other packages. Significant improvements in version 3 serve to mitigate. 
The integrated symbolic manipulation environment offered by DERIVE can outweigh its 
speed and programming disadvantage when non-numerical tasks are involved. 
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APPENDIX 1 

Utility file COUNT.MTH 

These functions may be used to produce a frequency count for elements in the vector v, 
given bucket boundaries in the vector b. 

INCV (v,p) increments the element of v. Next, we define functions to find the position of 
an element, e, of v. pos (e, v) produces an integer in the range [1,n], If e is less than the 
first element of v, pos (e, v) returns 1. If e is greater than the last element of v, pos (e, v) 
returns n+1. Otherwise POS (e, v) calls midpos (e, v) which returns an integer in the range 
[2,n], which indicates the bucket to which it belongs. 

POSITIONS (v, b) analyses each element in the vector v in this way. 

iNiT_C0UNT (v) initialises the count for each bucket to zero and count (v, b) maintains a 
tally of elements in vector v based on the class boundaries vector, b. 

INCV(v,p) :=REPLACE_ELEMENT (ELEMENT (v,p) +1, v,p) 

MIDPOS(e,v) :=MAX (ITERATES (IF (ELEMENT (v, i_) >=e, i_, i_+l) , i_, 1 , DIMENSION (v) ) ) 

POS(e,v) :=IF (e<ELEMENT (v, I) , I , IF (e>ELEMENT (v, DIMENSION (v) ) , 

DIMENSION(v) +1, MIDPOS (e,v) ) ) 

POSITIONS (v,b) : =VECTOR (POS (ELEMENT (v,u_) ,b) , u_, I , DIMENSION (v) ) 

INIT_COUNT (v) :=ITERATE (APPEND (w_, [0] ) , w_, [] , DIMENSION (v) +1) 

COUNT(v,b) :=ELEMENT (ITERATE ( [ I+ELEMENT ( c_, I ) , INCV (ELEMENT ( c_, 2 ) , 

ELEMENT (v, ELEMENT (c , 1) ) ) ] , C , [I, INIT COUNT (b) ] , DIMENSION (v) ) , 2) 


APPENDIX 2 

Utility file BARCHART. MTH 

The midpoints, left hand- and right hand end points of the bars to be produced are calculated 
by the functions x_mid (b) , x_endi (b) and x_end 2 (b) . These are combined in the 
vector x_ORDS (b) . rect (a,b,h,w) builds a rectangle of height h, and width w with (a,b) 
as centre point of its base. The function VBARS(ct,b,w) produces a vector of rectangles, 
using data in a vector ct, which can then be plotted. It is necessary to set the following plot 
options to obtain reasonable looking bars: Plot Color Cycling off (Color Auto No in DfD) and 
Points Connected and Size Small (in the State submenu in DfD). 

X_MIDS (b) :=VECTOR( (b SUB (i+I)+b SUB i) /2, i, I, DIMENSION (b) -I) 

X_ENDI (b) := (-b SUB 2+3*b SUB I)/2 

X_END2 (b) := (-b SUB (DIMENSION (b) -I) +3*b SUB DIMENSION (b) ) /2 

X_ORDS (b) :=APPEND (APPEND ( [X_ENDI (b) ] , X_MIDS (b) ) , [X_END2 (b) ] ) 

RECT (a,b,h,w) :=[ [a-w/ 2 , b] , [a-w/ 2 , b+h] , [a+w/ 2 , b+h] , [a+w/2 , b] ] 

VBARS (ct,b, w) :=VECTOR(RECT ( (X_ORDS (b) ) SUB i,0,ct SUB i,w),i,I, 

DIMENSION (X ORDS(b))) 
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Darstellung von Polyedern - 
Platonische und archimedische Korper im Mittelpunkt 

Hubert Weller, Lahnau, Germany 

Die platonischen Korper - The Platonic Solids 


Zur Bestimmung der Koordinaten der Eckpunkte 
werden auch die Dualitaten der Korper (z.B. Wiirfel 
<==> Oktaeder) ausgenutzt. Interessant ist die Be- 
rechnung der Koordinaten der Eckpunkte von Iko- 
saeder und Dodekaeder. 

For computing the vertices we use the duaiity 
of the solids (eg Cube Octahedron). Of spe- 
cial interest is the calculation of the vertices of 
an icosahedron and a dodecahedron. 

Icosahedron: 

A := [1,-s, 0], B := [1,s,0] and C := [s,0,1] with 
^ = AC = ^ = 2s gives s = — - - . 



Der Ikosaeder laBt sich aus einem Wiirfel aufbauen,wobei wir das Koordinatensystem in die Mitte des 
Wiirfels legen. Mit dem Ansatz A [1, -s, 0], B [1, s, 0] und C [s, 0, 1] laBt sich die Zahl s 
berechnen, da alle Kanten die Lange 2s haben. 

IKO_DODE.MTH 


#1: 

The I( 

:osahe 

idron in 

the 

cube; the cube: 





#2: 

[wl := 

[1, - 

1, -1], 

w2 :: 

: [1, 1, -1], w3 

:= [- 

1 

1 1 

% 
1 \ 

1 1 

1 

-1, 

-1]] 

#3: 

[w5 := 

[1, - 

1, 1], 

vj 6 := 

[1, 1, 1], w7 := 

[-1, 

1, 1], wS ;= [-1, -1, 

1]] 


#4: 

cube 

= [wl, 

w2 , w3 

. w4, 

wl , w5 , w6 , w2 , 

w 6 1 

w7, w3, w7, wS, w4, 

wS, 

w5] 



#6: 

cube 

■ proj 




#7: 

The \ 

/ertices 

of the 

Ikosaeder - 

The vertices of the Icosahedron 

#S: 

[a := 

[1, -s, 

0], b : 

= [1, s, 0], 

c := [s, 0, 1], d := [0, 1, s], e := [0, 1, -s]] 


#9: 


and so on 
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#10: ! 

Calculation of s with |ac| 

= 2s 

#11: 

2 

(a - c).(a - c) = (2.s) 


#12: 

1 ^5 


#15: ( 

22 ' 2 

Dnly the positive solution 

make: 

#14: 

JS 1 


#15: 

2 2 

ikosa = [b, a, c, b, d, c, 

d. e 

#16: 

ikosa := [b, a, c, b, d, c, 

d. e 

#17: 

ikosa-proj 



1 ■ 
2 . 



Auch der Dodekaeder laBt sich aus einem Wiirfel aufbauen, indem auf jede Wurfelflache ein 
Walmdach mit der Hohe h aufgesetzt wird. 



The dodecahedron can be con- 
structed from a cube attaching 
an Italian roof with height h on 
each side face. Let the length of 
the cube edges a = 2 and the 
length of all the dodecahedron 's 
edges s. So you can easily find 
values for h and s. 


Die Kantenlange des Dodekaeders sei s. Die Kantenlange des Wiirfels sei a = 2. Damit werden 
zunachst die Eckpunkte des Dodekaeders definiert: 

A ;= [0 ,2 ,2] , B := [-2 ,2, 2] , P := [- 0.5(2 - s), 2 + h, 1], Q := [-1 - 0.5s, 2 + h, 1] 
und R := [-1, 1 + 0.5s, 2 + h] usw... 

Da alle Kanten des Fiinfecks die Lange s (auch AP) und alle Diagonalen die Lange 2 haben (auch 
AQ), konnen die Werte fiir s und h berechnet werden. 
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IKO DODE.MTH (cont.) 


#1S: Der Dodekaeder urn den Wurfe! - The Dodecahedron around the cube 

#19: Das Walmdachdodekaeder - The hip- roof-dodecahedron 

#20: [cl ;= [0. 2. 2], c2 := [-2. 2. 2]] 

#21: [p := [- 0.5-C2 - s) , 2 + h, 1], q := [-1 - 0.5-s, 2 + h, 1], r := [-1, 1 + 0.5-s, 2 + h]] 

#22: and so on 

#23: Calculation of s and h; |clp| = s and |clq| = 2 


#24: (cl - p) ■ (cl - p) = s 


2 s - 4-s + S 2 
h + = s 


#23: (cl - q).(cl - q) = 2 = 


2 s + 4-s + S 
h + = 4 


#26: SOLUTIONS 


#27: SOLUTIONS 


^ 2 

2 s - 4-s + S 2 
h + = s , h 

V 4 

^ 2 

2 s + 4-s + S 
h + = 4, h 

V 4 



2 

7(3-5 + 4-s - S) 

2 

7(3-5 + 4-s - S) 

) 

1 

2 

2 



2 

7(- s - 4-s + S) 

2 

7(- s - 4-s + S) 

J 

1 

2 

2 


#2S: Only the positive solutions make sense! 


#29: SOLUTIONS 


2 2 
7(3.5 + 4.S - S) s - 4.Cs - 2)) 


[75 - 1 , - 75 - 1 ] 


#30: s ;= 7S - 1 


#31: h ;z 


7(_ 5 - 4.(5 - 2)) 


#32: h ;= 


4^5 


2 2 

#33: Eine Dodekaederf 1ache : 

#34: dof1l := [cl, p, q, c2, r, cl] 

#33: dofll-proj 

#36: [c3 := [-2. 0. 2], c4 := [0. 0. 2], c3 := [0. 2. 0]] 

#37: [c6 := [-2, 2, 0], c7 := [-2, 0, 0], cS := [0, 0, 0]] 

#3S: cube2 := [cl, c2, c3, c4, cl, c3 , c6, c2, c6, c7, c3, c7, cS, c4, cS, c3] 


cube2 ■ proi 


Next page shows the plot of the cube together with one face of the dodecahedron. 
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Verdeckte Kanten und Determinanten Hidden lines and Determinants 

Die bisher benutzte Darstellung ist besonders einfach, denn ein Korper wird durch eine einzige Punkt- 
folge (d.h. durch eine Matrix) beschrieben. Allerdings hat diese Form der Darstellung auch groBe 
Nachteile, denn u.U. werden einige Kanten mehrmals gezeichnet und man sieht alle verdeckten Kan- 
ten. Deshalb haben wir uns iiberlegt, wie die Form der Darstellung verbessert werden kann. 


Wenn sich ein konvexer Korper im Raum 
bewegt, weiB man, dass eine Flache, die vorher 
sichtbar ist, in dem Moment, in dem sie 
„verschwindet“ (oder eine vorher nicht 
sichtbare Flache „auftaucht“), ihre 
ORIENTIERUNG andert. Dies ist die ent- 
scheidende Idee bei den folgenden Uberlegun- 
gdn. DETERMINANTE einer 2x2 -Matrix laBt 
sich deuten als Flacheninhalt des von den 
beiden Zeilen-vektoren aufgespannten 
Parallelogramms. Das Vorzeichen gibt dabei 
eine Information iiber die Orientierung, d.h. 
wenn der erste Zeilenvektor auf direktem Weg 
in den zweiten gedreht wird, dreht man bei 
einer positiven Determinante im mathematisch 
positiven Sinn, bei einer negativen 
Determinante im mathematisch negativen 
ZifiQ.ist 


Das von den Vektoren a = [3 l] und 
^ = [l 2] aufgespannte Parallelogramm hat 
den Flacheninhalt 5. 


It is known that a side face of a convex 
solid when moved in space is changing its 
orientation at the moment when it is 
changing its visibility. This is leading to the 
following algorithm: 

The determinant of a 2 x 2 matrix can be 
interpreted as the area of a parallelogram 
built by the two row vectors r^ and r 2 . Its 
sign gives information about its orientation. 
(= 1 if a positive rotation would move r^ 
into r 2 ). 

The area of the parallelogram built by the 
vectors a = [3 l] and b = [l 2] is 5. 
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Das von den Vektoren a = [-l 3] und b =\2 2] 

aufgespannte Parallelogramm hat den Flacheninhalt -8; 
dieser ist negativ, weil das Parallelogramm eine andere 
Orientierung hat. Wir benutzen das als 
Entscheidungskriterium dafiir, oh eine Flache 
gezeichnet werden soil oder nicht. 

Here the parallelogram's area is -8. We use the 
sign as a criterion whether a side face should be 
plotted or not. 



Eine andere Art der Darstellung von konvexen Korpern 

Another way to represent convex solids 

Die Eckpunkte werden wie bisher als Zeilenvektoren beschrieben. Z.B. bei einem Wiirfel: 
pl: = [l,0,0], p2: = [l,l 0], 

Da allerdings bei jeder Flache entschieden werden soil, ob sie gezeichnet werden soil oder nicht, 
stellen wir eine Flache als Punktfolge (geschlossener Streckenzug) dar: 
f 11 : = [pi , p2 , p6 , p5 , pi ], f 12 : = [p2p3p7p6p2 ], ... 

Bei der Definition der Flachen muB darauf geachtet werden, dass die Punktfolge beim Blick von 
auBen auf die Flache im mathematisch positiven Sinn durchlaufen wird. 

Einen Korper beschreiben wir nun als Vektor von Flachen, d.h. einen Vektor von Matrizen: 
cube:=[fll, fl2, fl3, fl4, fl5, fl6] 

Zum Test, ob eine Flache gezeichnet werden soli oder nicht, werden nur die ersten 3 Punkte dieser 
Flache benutzt. Das Vorzeichen der Determinante derjenigen Matrix, in deren ersten Zeile die Pro- 
jektion des Vektors vom 2. zum 3.Flachenpunkt und in deren zweiten Zeile die Projektion vom 2. zum 
1 .Flachenpunkt steht, liefert das Kriterium, ob die Flache gezeichnet wird. (Den i-ten Punkt der k-ten 
Flache eines Korpers erhalt man mit dem Befehl cube SUB k SUB i oder iiber 
ELEMENT (cube, k, i) ) . 

The points are written as row vectors. Each side face is defined by a closed sequence of 
points. The solid is described by a vector of side faces (see cube : = [ . . . ] above). 

Defining the side faces (polygons) you have to take care that the sequence of points has to 
run through in positive direction when you look at the polygon from outside of the soiid. For 
testing the orientation you need oniy the coordinates of the first three - projected - points. 


TESTCsolid, 


k) ;= DETC 


solid 


-solid , solid 


- solid 


■proj) 


L k,3 k,2 k,l k,2j 

Eine Flache wird gezeichnet, weim beim Test ein positiver Wert berechnet wird, sonst wird sie nicht 
gezeichnet! 

The side face wili be piotted if the test results in a positive vaiue! 


PIC_2_D(solid) := VECTORClFCTESTCsolid, k) > 0, solid -proj), k, 1, 

k 


DIMCsol id)) 
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HID_LINE1.MTH 


#1: Verdeckte Kanten bei konvej<en Kdrpern - Hidden lines on convex solids 

#2: Die Determinante a1s orientierter Flacheninha1t wind benutzt 

Die Eckpunkte des Kdrpers - the vertices 
#4: [pi := [1, 0, 0], p2 := [1, 1, 0], p3 := [0, 1, 0], p4 := [0, 0, 0]] 

#5: [p5 - [1. 0. 1], p6 - [1. 1. 1], p7 [0. 1. 1], pS [0. 0. 1]] 

#6: Die Flachen - dabei auf mathematisch positiven Umlaufsinn achten ! ! 

#7: The faces - take care of the positive rotation direction 


#S: [fll 

:= [pi, 

P2, 

p6, 

P5, 

pl] , 

fl2 ;= [p2, 

P3, 

P7. 

p6, 

p2] , 

fl3 ;= [p3, 

P4, 

pS, 

P2, 

1 1 

1 1 

nn 

Cl 

1 1 

-h 

:= [p4, 

Pl, 

P5, 

pS, 

p4], 

fl5 ;= [pS, 

P5, 

p6, 

P7, 

pS], 

fl6 := [p4. 

P3, 

P2, 

pl, 

p4]] 


#10: Der Kdrper - the solid 

#11: cube := [f1l, fl2, fl3. fl4, fl5. fl6] 

#12: Die Projektionsmat r ix - the projection matrix 


#13: proj 


-0.5 -0.25 ■ 

1 0 

0 1 


#14: Der Test, ob eine Flache gezeichnet werden so11 - plotting test: plot or not plot? 

TEST (solid, k) := DETCFsolid -solid , solid -solid 1-proj) 

#15: L k,3 k,2 k,l k,2j 

#16: Das 2C^Bi1d - the 2C^projection 

PIC_2_D(so1id) := VECT0RClF(TESTCso1id, k) > 0, solid -proj), k, 1, DIH(solid)) 

#17: k 

#18: IdiiMiMaBHl 

Plot #18 in order to obtain an image of the cube: 



Bewegungen des Korpers im Raum konnen mit 
Hilfe von Matrizen dargestellt werden, dabei 
werden ausschlieBlich die sichtbaren Flachen 
gezeichnet. Als Beispiel soil ein Wiirfel auf die 
Spitze gestellt werden. 


Movements in space can be realized using 
matrices. Only the visible side faces will be 
plotted. As an example we represent a 
cube standing on its top. 
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#19: Wir stellen einen Wurfe! auf die Spitze! - The cube on its top! 

#20: 45°-Drehung urn die z-Achse - 45°-Rotation around z-axis 

COSCffl) SINC«) 0 


#21: DZ(ffl) ;= 


- SINCh) COSCh) 0 
0 0 1 


cubel := VECTORCcube .02(45.1°), k, 1, DIH(cube)) 

#22: k 

#23: Rotation around the x-axis - Drehung urn die x-Achse 

#24: The cube's diagonal is rotated into the z-axis (around the x-axis) 

10 0 


#25: DX(ffl) := 


0 COS (a) SIN(ffl) 

L 0 - SIN(ffl) COS(ffl) J 

^ 1 > 


#26: p := ATAN 






r TT '] 


#27: cube2 VECTOR 

cubel .OX 

P 

, k, 1, DIM(cubel) 


L k 

L 2 J 

J 


#2S: Der Wiirfel steht am Kopf! - The cube is standing on its top! 


#29: axl 


" 1 

0 

0 ■ 


■ 0 

0 

0 " 


■ 0 

0 

0 " 

2 

0 

0 


2 

0 

0 


2 

0 

0 

Sl 

00 

00 


00 

00 

00 


00 

00 

00 

0 

1 

0 

ax2 :z 

0 


0 

;= 

0 

0 

0 

0 

2 

0 


0 

2 

0 


0 

2 

0 

00 

00 

00 


00 

iS> 

00 


00 

00 

00 

0 

0 

1 


0 

0 

0 


0 

0 


_ 0 

0 

3 _ 


_ 0 

0 

3 _ 


_ 0 

0 

3 . 


#30: [axl.proj, PIC_2_D(cube) ] #32: [ax2.proj, PIC_2_D(cubel) ] #34: [ax3-proj, PIC_2_D(cube2) ] 



p 36 


Hubert Weller: Regular Polyhedrons - Hidden Lines 


D-N-L#29 


Literatur - References 

H.-J. Elschenbroich/J.-C. Meiners, Computergraphik und Daestellende Geometric im Unterricht der 
Linearen Algebra, Diimmler (Bonn) 1994 

W.Kroll (Hrsg.), Neue Impulse fiir die Raumgeometrie, mathematik lehren Heft 77 1996 
Friedrich Verlag (Velber) 

H. Weller, Determinanten in einem Kurs Lineare Algebra, DdM 1 , 1979 , (S. 62-72) 

E.-C.Wittmann, Elementargeometrie und Wirklichkeit, Vieweg (Braunschweig) 1987 



In DNL#21 Richard Schorn presented the "Bucky Ball" - a model of the C60-molecule (12 pentagons 
and 20 hexagons. I used the data of C60LAPHI.MTH and Richard Schorn 's sketch from page 42 to 
present the Bucky Ball and let it rotate round the z-axis. Josef (BUCKY.MTH) 
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Dear Joseph, North Balwyn, 5 April 1996 

I have various findings, which I believe both the Derive community and the mathematical community will find of 
interest, so I shall start to send them to you successively, this letter being the largest by way of introduction. 


The first, I have called the "Cesaro Glove-Osculant", a Super-Osculant, eine besondere Schmiegungstan- 
gentenkurve, a sheath osculant. It opens a new and innovative look at osculants, and simultaneously gives one the 
power to generate images that closely simulate a sought-after image and that may have been too difficult to discover 
by other means. It was introduced by Ernesto Cesaro in his famous tome, Vorlesungen Uber Naturliche Geome- 
trie around the turn of the century. 

I have opened a mathematical 'can of worms', too much for one person to develop to its full potential, so I invite 
fellow members of DUG or readers of the Newsletter to write to me, in English, and perhaps we can form a team for 
further development of this line of research. 

Enclosed the paper SUPEDUPE.WPF on disk and print-out, as well as Derive files, OSCCESln.MTH and two image 
files, TAUTUREn.WPG on disk. 

Yours sincerely, 

David Halprin. 


SUPER DUPER OSCULANTS 

David Halprin, North Balwyn, Australia 

ABSTRACT & INTRODUCTION 

As an initial analogy, consider the case of a form-fitting glove (Handschuh), [or 
sheath (Scheide)], for a hand (or surface) in three dimensions. Suppose that there 
exists a machine for manufacturing the sheath, which depends on a set of three 
parametric equations being fed into the computerised instructions and, if one 
provides these equations to define the surface, then out comes the sheath exactly 
form-fitting. Suppose there exists a set of surfaces, whose definition is in terms of 
intrinsic variables, and which cannot be transformed into a Cartesian parametric set. 
This means that the machine cannot produce the appropriate sheaths until and 
unless an algorithmic method of close approximation can be devised to transform the 
equations appropriately. 

This paper deals with the two-dimensional counterpart to the sheath, which we 
can describe as a piece of bent wire, which can be laid on top of the curve in the 
plane, and covering it exactly. If the wire were bent by a machine, which depended 
on the equation of the curve being given in a pair of parametric Cartesian equations, 
then it would produce an exact replica of the curve. If, however, we could only 
represent the curve in Intrinsic variables, then we would hope to approximate the 
curve with a pair of parametric equations, which, although not being exactly the 
defined curve, would so closely approximate it as to be almost error free for most of its 
chosen arc. Such a curve could be described best as a very special type of osculant, 
as yet, not dealt with in papers and texts on osculants, since it osculates over almost 
all the curve, rather than some small arc centred about a nominated point, or group 
of points, on the curve. 

Hence this paper is to show the theory and derivation of the algorithm and 
then to exemplify it in a routine, (in Derive or DeriveXM), which can illustrate almost 
any curve one desires to show. (See enclosed OSCCESIN.MTH files) 
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SUPER-OSCULANTS (THEIR KITH AND KIN) 

An Osculant in the plane can be looked upon as one of several possible 
mathematical devices or tools, (some inter-related), to assist one in the investigation 
of a geometrical image of all, (or part), of a nominated base curve in the plane. Its 
usage may be:- 

A) direct, in any one of at least four separate categories of osculants, 
(Penosculants, General Osculants, Chief-Osculants & Super-Osculants) 


or 

B) indirect and/or incipient, (perhaps via overt consideration of order-of-contact), 
and/or one or more of the ascending series of geometrical properties, 
(qualitative and quantitative), starting with:- 

1) Gradient, [first order], then 

2) Radius of curvature, [second order], and 

3) Deviation (aberration) or spiralation, [third order] . 

4) ^ not to mention often limitless higher order possibilities. 

A plane curve may be investigated in many different ways, to elucidate its 
various geometric properties. One way draws comparison of the base curve with 
other curves, (such as osculants), by comparing features, viz:- 

FIRST ORDER 

The gradient of the base curve at a nominated point is the gradient of the 
osculating straight line, (the tangent), at that point. 

SECOND ORDER 

The radius of curvature of the base curve at a nominated point is the radius of 
the osculating circle. 

THIRD ORDER 

The deviation (aberration) of a base curve at a nominated point is the deviation 
of an osculating parabola. 

The spiralation of a base curve at a nominated point is the spiralation of an 
osculating logarithmic spiral. 

In all, (but one), of the above-cited cases we have used, as nominated 
osculants, those curves, (called "Chief-Osculants"), which have the chosen quality as 
a constant, namely: 

a straight line has a constant gradient, 
a circle has a constant radius, and 

a logarithmic spiral has both constant deviation and constant spiralation. 

N.B. A parabola does NOT have constant deviation, despite some authors having 
errone-ously stated an alleged constancy. 
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The chief-osculants of higher orders than 3 are difficult to plot, especially from 
the fifth order onwards, hence the particular geometric image associated with each is 
nigh impossible to envisage, yet a mathematician can deem it to be so visualised for 
the purpose of using it as a referential device, (both qualitatively and quantitatively), 
when describing some new base curve. There are, however, better ways to gain 
advantage from the osculant concept, namely by using: - 

A) Some specially selected general osculants of a chosen order of contact, 
[requiring integration, (either producing a closed form resulting equation, or an 
approximated plot by numerical integration techniques)], or 

B) Some specially devised super-osculant, which requires no integration at all, but 
will require differentiation of the intrinsic equation of the base curve up to the 
nth. order to produce an osculant with an (n-l)th. order-of-contact. This is 
often so much like the base curve, between the selected plotting limits, that, to 
all intents and purposes, it may well be looked upon as the base curve's image, 
morphologically speaking, and thereby one has obtained a geometric image of a 
curve, that otherwise may be almost impossible to plot or at least most difficult 
to plot accurately. 

Now to demonstrate the techniques for A) and B) above. 

A) There is a recursive relationship that explicitly defines all such Chief-Osculants. 

RECURSIVE RELATIONSHIP 

With respect to those Intrinsic Qualities of the Chief-Osculants, numbered 
above with the numerals only, in brackets:- 
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The establishment of these Chief-Osculants in differential form, (as illustrated 
in the Table at end of this paper), does not necessarily mean that we can EASILY 
identify either the Curves or the Geometric Images, so portrayed, since the Calculus 
involved for some of these probably involves indefinite integrals, which cannot be 
represented in terms of known functions, but can be drawn by evaluating by means 
of numerical approximation, or on a computer monitor, with a 'graph-plotter' 
program, such as Derive or DeriveXM, which integrates as it plots each point. 

B) The Super-Osculants are derived from one or more of the six Intrinsic 
Equations to a particular base curve. Even though there are six possible pairs of 
parametric equations for such a Super-Osculant, the ultimate choice, that produces 
the greatest semblance to a major section of the base curve, is most likely to come 
from either of the Whewell Equations, since, of the three intrinsic variables, [the 
radius of curvature (r rho), the arclength (s) and the tangential angle (f phi)], r being 
in an equation, (as it is in the two Cesaro equations and two Euler equations), is 
better avoided where possible, especially if it is the independent variable, since, either 

1) as a parameter in a pair of equations, it cannot be expected to range between two 
nominated extremes as the other two variables do so readily, or 

2) it often occurs in the denominator of the terms of the Taylor series and it 
necessarily takes a zero value at cusps and poles. 

This is an overview of what Cesaro did to obtain the equations of his osculants, 
and then my own extensions of his ideas, which, perhaps, he did not consider, 
opining, perhaps, that what he did was all there was to do. 

He assumed that there was no ambiguity, let alone more than two possible ways 
of defining order-of-contact. This assumption of his was implicit, since he never 
stated an 3 dhing on the matter, other that attest that the order of contact was defined 
by the evolutes of increasing order. In other words, he states that with a second 
order osculant, the osculation was representative of a commonality between base 
curve and osculant, that they shared the same first order contact, namely gradient. 


Now will follow a couple of pages, some of them full with large mathematical expressions, an 
interesting table and two graphs which should illustrate the power of the Super Osculants. I 
hesitated to publish this paper because it looked too "heavy" and I wrote back to David asking 
for some additional explanations and more illustrations. No answer. But in this January together 
with his membership renewal arrived a letter from Australia with David's complain, why his 
paper has not been published until now despite the announcements on all pages #2 in many 
DNLs. I answered and asked once more for some additional comments. Now I received the 
answer very soon and I find it worthwhile to be published below. It will support you 
understanding my requests and his ideas as well. 






D-N-L#29 

David Halprin: SUPER DUPER OSCULANTS 

p 41 



Dear Joseph, North Balwyn, 10 February 1998 


I received your letter today; many thanks for your prompt reply. Unfortunately I have not received your earlier 
letter, requesting more graphics, so I am rectifying that immediately. I was foolish to omit them since the main point 
of the Glove-Osculant algorithm is to demonstrate how good it is for approximating major arcs, (at least), of many 
common curves, as well as many curves, which have no simple closed-form algebraic expression. I have included 1 1 
such images on the enclosed disc as well as an up-dated version of the file that achieves that end, OSCUWHEl.MTH 
as well as the OSCUWHEl.PRT file for your easy reading. I have inserted relevant comments about the images in 
those files, and I shall print them out, together with all 1 1 images in the accompanying letter DUGIMAGE.WPF. 

Because you told me that you found nothing on Cesaro to date, I include some eulogy-type articles in Italian, 
written by people, who knew and respected the man and his researches into many disparate, (and sometimes 
untouched), areas of mathematics. Unfortunately I do not speak Itahan, therefore I could only glean a smidgen of 
their contents. Hopefully you will gain more than me from the papers. 

Your query re so many equations on pages 6-8 is very pertinent, especially if you have to retype them. I hope 
you have access to the WordPerfect Word Processor V 5.1 or later, since all those equations are created by the 
equation editor. In summary, I put in all the details so that the incredulous reader(s) won't have to take my word for 
it, perhaps thinking that I am using a vivid imagination, rather than an in-depth analysis of Cesaro's work with further 
developments by myself. Whereas he only considered this method for those equations, which I labelled Cesaro 
Type-1 p = C(s) there are 5 other forms of Intrinsic Equations, and each one had to be tackled and solved by me 

taking many tens of hours, by hand, before I eventually was able to develop an algorithmic routine using declared 
functions in Derive and thereby enabhng not only a valuable time-saving approach, but also no apparent limit to 
which high order-of-contact the glove-osculants make with the base curve. I sometimes went to the 15th order, which 
imphed 15 successive differentiations of these declared functions in a canonical form, so that when the result was 
obtained, it could be used for any nominated base curve by slotting it into the general form by using the Derive 
Simplify Command. That is why the .MTH file has thousands of terms in the vector [X(s),Y(s)]. 

I shall also include on the disc the two .WPF files, being this letter DUGllDLH.WPF and the accompanying 
letter DUGIMAGE.WPF 

I take the hberty of quoting from my mail of 1996, wherein I request some readers to contact me, hopefully :- 

/ have opened a mathematical 'can of worms', too much for one person to develop to its full 
potential, so I invite fellow members of DUG or readers of the Newsletter to write to me, in English, and 
perhaps we can form a team for further development of this line of research. 

Yours sincerely, 

David Halprin. 


To give an impression of David's contribution I add two screen shots of super osculants: 



The NEPHROID and its osculant 


The ASTROID and its osculant 


p 42 


Reactions on TITBITS 


D-N-L#29 










D-N-L#29 


Reactions on TITBITS 


p 43 


2 Implementation 

A single line in DERIVE based on (2) defines a function which is equivalent to J.Wiesenbauers's function 

GEN_LUCAS(...): 

FPL.MTH 

#1 : R (gv_, n_, sv_) —ELEMENT ([[0,1], gv_] ^ (n_-l) ^sv_, 1) 

#2 : FIBONACCI (n_) : =R ( [ 1 , 1 ] , n_, [0,1]) 

#3 : VECTOR (FIBONACCI (k),k,l,10)=[0,l,l,2,3,5,8,13,21,34] 

#4: PELL(n_) : =R ( [ 1 , 2 ] , n_, [0,1] ) 

#5 : VECTOR (PELL (k),k,l,10)=[0,l,2,5,12,29,70,169,408,985] 

# 6 : LUCAS (n_) : =R ( [ 1 , 1 ] , n_, [1,3] ) 

#7 : VECTOR (LUCAS (k),k,l,10)=[l,3,4,7,ll,18,29,47,76,123] 

#1 looks considerably simpler as J.Wiesenbauer's definition of GEN_LUCAS(...) in [1]. The price for this 
simplicity is speed. The following table compares computational times for the functions in NUMBER.MTH and 
the functions #2, #4 and #6 based on the matrix method: 


Function 

NUMBER.MTH 

Matrix method 

M.m with DfD 4 

PELL( 100000) 

149.5 sec 

12.2 sec 

6.9 sec 

LUCAS(IOOOOO) 

1.6 sec 

3.6 sec 

1.9 sec 

FIBONACCI(IOOOOO) 

2.2 sec 

3.6 sec 

1.9 sec 


The computation times refer to a DfW 4.0 and a PENTIUM 133 with 32MB RAM. (I added the DfD times, J.) 

EIBONACCI(...) is defined with two calls of LUCAS(...) in NUMBER.MTH, and LUCAS(...) is basically 
the same function as GEN_LUCAS(...) in [1]. This explains the rise from 1.6 to 2.2 seconds in the second 
column. There is no difference in the matrix version, because only the start vector is different, but the matrix is 
the same. The NUMBER.MTH version of PELL(...) uses recursion, and this fact explains the 149.5 seconds in 
contrast to the 12.2 seconds of the matrix method. 

3 Conclusions 

1. The simple and transparent matrix method is roughly half as fast as Wiesenbauer's method, but 
considerably faster than the recursion. The matrix method works even in the case where the discriminant 
of the characteristical polynomial is zero. 

2. The matrix method allows generalizations to recursive sequences of higher order. Eor example: Let {a„} 

be defined by a\ = x\, = X 4 and = a + P + y .. 2 + 6 1 then 


^ 0 

1 

0 

0" 

n-\ 

4 V 3 

Xi 

( a \ 

0 

0 

1 

0 


•^2 

««+i 

0 

0 

0 

1 


X3 

«« + 2 

V a 

p 

y 

6y 


7 xj 

V + 3 J 


It is evident how to generalize further. 

3. It is possible to construct formulas for direct computation of the recursive sequences under consideration, 
which generalize the Binet-formula for the Eibonacci sequence. The generalized Binet-formulas are 
appropriate linear combinations of n-th powers of eigenvalues of the corresponding matrices [2] 

References 
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[2] Welke, S.: Line Verallgemeinerung der Binet-Formel fiir linear rekursive Folgen, to appear in MNU 
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/ informed Johann about this paper and he asked me to fax it immediately to his institute (Vienna 
Technical University). He e-mailed an answer containing some comments and sent a Cc to Stefan 
Welke. Stefan expressed his thanks, considered his advice and made some changes as: 


The table should read: 


The following table shows some examples: 


Sequence 

a 

b 

a 

p 

Fibonacci 

1 

1 

1 

1 

Pell 

0 

1 

1 

2 

Lucas 

1 

3 

1 

1 


That is the improved FPL.MTH 

#1: R (gv, n, sv) = ( [ [ 0 , 1 ] , gv] ^ (n-1 ) ’^sv) SUB 1 

#2 : FIBONACCI (n) :=R([l,l],n, [1,1]) 

#3 VECTOR (FIBONACCI (k),k,l,10)=[l,l,2,3,5,8,13,21,34,55] 

#4: PELL(n) :=R( [1,2] [0,1] ) 

#5 : VECTOR (PELL (k),k,l,10)=[0,l,2,5,12,29,70,169,408,985] 

#6: LUCAS (n) :=R( [1,1] ,n, [1,3] ) 

#7 : VECTOR (LUCAS (k),k,l,10)=[l,3,4,7,ll,18,29,47,76,123] 

Stefan writes that he - in contrary to Johann - is finding the use of both auxiliary functions and long 
variable names as an advantage. He prefers the shifting of the index, because he intends to present a 
general solution. Johanns suggestion uses an intial value of index 0, which has to be known in each 
case. 

Find below Johann 's concluding message: 

Dear Josef, 

Thank you very much for sending me the interesting paper of Stefan Welke. Though due to my notorious 
idiosyncrasy to auxiliary functions I would have preferred the definitions 

FIBONACCI(n) := ([[0, 1], [1, l]]^n*[0, 1])SUB1 

LUCAS(n) ([[0, 1], [1, l]]"n*[2, 1])SUB1 

PELL(n) ([[0, 1], [1, 2]]"n*[l, 0])SUB1 

I don't object to his original routines either which are short and sweet as well and, above all, incredibly fast! As a 
matter of fact, they are more than ten times faster than former DERIVE-implementations based on the original 
recursion and only about two times slower than the current versions in DfW 4.08. Moreover Stefan Welke is 
right in saying that you don't have that restriction about the nonzero discriminant (even if this condition is 
usually fulfilled!) 

This raises the question why I didn't take that method into consideration. The answer is quite simple: In many 
important applications you need the values of those series mod N for huge indices (say with 100 digits; see the 
example in DNL #26, p27, to see what I mean!) and here the matrix-method simply can't keep up. (You have to 
develop a modulo-arithmetic for matrices to avoid a memory overflow and all the performance goes to blazes!) 

Anyway, concluding I would like to say that whenever didactic aspects are more important than absolute top 
performance Stefan Welke's routines are certainly an alternative one should consider. In particular, the people 
who wrote the chapter about recursive functions in the DERIVE-manual should study his paper very carefully. 

Cheers, Johann 
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Another letter arrived - "Bits and Pieces in Number Theory" - from Andreas Eder, Munchen, 
Germany: 


Dear Josef, 

first let me express my thanks for your good work. I just received DNL #28 and had a short look over it. It 
seems to be full of interesting articles once again! 

Now to my little contribution: I especially like the Titbits from Algebra and Number Theory' by Johann 
Wiesenbauer, partly because I myself am very interested in Number Theory, and because I constantly learn new 
programming tricks by trying to understand his programs. So I discovered what he calls the ITERATE( , , , 1) 
trick. A very clever thing, but to be honest, I don't think it makes the programs very readable. So I thought 
about it and tried to encapsulate that trick in a bit of 'syntactic sugar'. This is what I have come up with: 


LET(v,e,b):=ITERATE(b,v,e,l) 


which should be read as: let the vars be bound to expressions in body (like the let construct in Lisp or Scheme). 
I think it makes programs more succinct and readable, e.g. 


LIN_DIOPHANT(a,b,c):=LET([d_,uJ,GCD_EX(a,b),IF(MOD(c,d_)=0,u_*c/d_+[b*@,-a*@])) 


This version of LIN DIOPHANT also introduces my version of EXTENDED GCD, which I called GCD EX. 


GCD_EX(a,b):=LET([x_,g_,y_,m_],ITERATE(IF(MOD(u3_, v3_)=0,[u3_,v3_,ul_,vl_],^ 

[v3_, MOD(u3_,v3_),vl_, ul_-FLOOR(u3_, v3_)*vlj),[u3_,v3_,ul_,vl_],-^ 

[a,b,l,0]),[g_,[m_, (g_-m_*a)/b]]) 

You can see threee differences with respect to EXTENDED_GCD: first it uses LET, second it uses named 
variables instead of a vector and subscripts and third and most important - it does away with the components 5 
and 6 of the vector in EXTENDED_GCD. These are not really needed, because if g = m*a+p*b, knowledge of g 
and m is sufficient to calculate p. This makes the program both shorter and more readable, and it should be 
slightly faster as well. 

In all these, and some more cases my LET construct worked as expected. I hope it will be of some use to others, 
as well as my streamlined version of EXTENDED_GCD. 

Yours, Andreas Eder Andreas.Eder@t-online.de 

From Jan Vermeylen, Kapellen, Belgium: 

Today I received DNL#28. Couldn't stop reading it! Very interesting articles indeed, (as always, but still ....). In 
the 4.07 version of the NUMBER.MTH utility file you can find a new function SORT(a) which sorts the 
elements of a vector with different rational numbers. The function that I incorporated in my STAT.MTH file can 
also sort vectors with non-different irrational numbers. Check the attached file: I made two versions SORTX and 
SORTY (with and without frequency) - Best greetings, Jan. 

(You will find both functions SORTX and SORTY in the file RESORTMTH, Josef) 

Do spread the news registered DfW users can have free upgrades from the SWH site now! 

The TITBITS master's Johann Wiesenbauer's answer: 

Sorry for the delayed answer caused by my holidays. I assure you that I'm as pleased as Punch about every 
feedback I get from readers of my 'Titbits'. As an avid programmer I can vividly imagine that you were very 
enthusiastic about your beautiful SORT-routine when it finally yielded the correct results. Therefore I hate to 
inform you that it has a serious flaw: Its asymptotic growth is like 0(n'^2), where n is the length of the list in 
question, which is far too slow for practical purposes. (As a matter of fact, for lists with more than 20 elements 
your SORT-routine becomes even slower than that of Sergey Biryukov speeded up by me to no avail in DNL 
#23, p39.) In other words, if you still feel like sorting with DERIVE, you are strongly advised to change to a 
better sorting algorithm like e.g. quicksort (here the corresponding magnitude is 0(n log n) on average, i.e. its 
growth is almost linear!) 

In the meantime, I can offer you a general SORT-routine of my own that should work in most cases. (It is longer 
than necessary due to two tiny DERIVE-bugs which I have described in the attached file. I am sending a carbon 
copy of this e-mail to A1 Rich hoping that he has mercy on us.) 

Please, keep me informed about your future programming efforts related to the topic above! - Cheers, Johann 
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RESORT.MTH 

SORTX (v) := (ITERATE (IF (DIMENSION (x SUB I ) =0 , x, [ SELECT ( k/=MIN (x SUB l),k, 
X SUB I), APPEND (x SUB 2,[MIN(x SUB 1 )])]), x, [v, []]) ) SUB 2 

SORTY(v) :=APPEND (VECTOR (VECTOR(x,i, I, DIMENSION(SELECT(k=x,k,v) ) ) , x, 
SORTX (v) ) ) 

"(Your sort-routines)" 

"Comparing some DERIVE-routines for sorting of numbers" 

" ( (c) Johann Wiesenbauer, Vienna) " 

SORTO (v) :=LIM(TERMS (v*VECTOR (x_^ (-k_) , k_, v) ) ,x_, 1) 


(SORTO(v) is very fast, but works only if v is a list consisting of different(!) nonzero rational numbers.) 

SORT(v) :=ITERATE (IF (v_ SUB DIMENSION (v_) =0 , APPEND ( SELECT (u_<0 , u_, v_) , 

SELECT (u_=0,u_,v_) , SELECT (u_>0 , u_, v_) ) , v_) ,v_, VECTOR (v SUB k_, k_, 
LIM (DIF (TERMS (EXPAND (VECTOR (x_^APPROX (-v_) , v_, v) *VECTOR (y_^k_, 
k_,l, DIMENSION (V) ) ) ) ,y_) , [x_,y_] , [1,1] ) ) ,1) 


You should use sort(v) if the preconditions above are not fulfilled. 

SSORT_AUX(v,n,p) :=VECTOR(IF(MOD(i+p,2)=0,IF(i=l OR v SUB (i-l)<=v SUB i, 

V SUB i,v SUB (i-1) ) , IF (i=n OR v SUB i<=v SUB (i+1), 

V SUB i,v SUB (i+I))),i,n) 

SSORT (v) :=ITERATE ( SSORT_AUX ( SSORT_AUX (w, DIMENSION (v) , 0) , DIMENSION (v) , 1) , 
w, v) 


The sort-routine by S.Birynkov speeded np by me in an earlier DNL: 

[a:=[l, 1,2, 1,2, 3, 1,2, 3, 4, 1,2, 3, 4, 5], b:=[pi, SQRT (2) , LOG (2) , 0, LOG (2) ] ] 

SORT (a)=[l, 1,1, 1,1, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5] 

SSORT (a)=[l, 1,1, 1,1, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5] 

SORTX (a) =[1,2, 3, 4, 5] 

SORTY (a)=[l, 1,1, 1,1, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5] 

SORT (b) = [0, LN (2) ,LN(2) , SQRT (2) ,pi] 

SSORT (b) = [0, LN (2) ,LN(2) , SQRT (2) ,pi] 

SORTX (b) = [0, LN (2) ,SQRT(2) ,pi] 

SORTY (b) = [LN (2) ,SQRT(2) ,pi] 

SORTY (b) = [0, LN (2) ,LN(2) , SQRT (2) ,pi] 

There has been no difference so far, right? 

Bnt look at the following important example from nnmber theory: 

FAREYO(n) :=( ITERATE ( IF (v_ SUB l=v_ SUB 2 AND v_ SUB 2=n,v_, [IF(v_ SUB 1= 

v_ SUB 2,l,v_ SUB l+l),IF(v_ SUB l=v_ SUB 2,v_ SUB 2+1, v_ SUB 2), 
IF(GCD(v_ SUB l,v_ SUB 2 ) =1 , APPEND (v_ SUB 3, [v_ SUB l/v_ SUB 2]), 
V SUB 3 ) ] ) , V , [1,1, [ ] ] ) ) SUB 3 


You should simplify each of the following lines (inclnding the next one!!) and compare! 


[ c : ==FAREY0 ( 2 0 ) , d : ==FAREY0 (40)] 


[DIMENSION (c) , DIMENSION (d) ] 
SORTO (C) 

SORT (c) 

SSORT (c) 
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SORTX (c) 

SORTY (c) 

SORTO (d) 

SORT (d) 

SSORT (d) 

SORTX (d) 

SORTY (d) 

My times (in s on a Pentium 200 PC): 0.1,0.4,17.3,14.9,17.3 and 0.6,1.6 


You might remember "Problem 4 of the CAS Competition (DNL#27, p25). Johann tried an 
approach to that challenge using the TI-92: 

A TI-92 Approach to Problem 4 

The following TI-92 function returns the n-th coefficient of the polynomial in question scaled 
down by the factor E-500. Note that the calculator must be in approximate mode! 

coef f ( n) 

Fu nc 

Local a , b, i 

{1. e- 500, 0,0, 0,0, 0,0}^ a 
- { 3 5 0 0, 1 0 0 0 5, 1 6 5 1 2, 1 9 5 1 8, 2 0 0 2 0, 1 4 5 
For i , 1 , n 

left(augment((dotP(a,b)/i},a) 
b + {3, 5, 6, 6, 5, 3, l}^b 
EndFor 
Return a [ 1 ] 

EndFunc 


Computing coeff(3000) yields 
3.9739422655873e926 after about 5 min, i.e. 11 
significant digits of the number in question! 
(Note that 500 must be added to the exponential 
part.) 


/ tried the "official" solution presented on TTs home page prob4(), consumed one set 
of batteries the first night running the program and started once more with a fresh set, but I 
interrupted after 13 hours. The I/O screen showed some progress: 12, 24, , 300, 600. 

In the directory I found three lists for intermediate results containing approximately 16000, 
45000 and 47000 bytes. Josef 


1 5 , 6 0 0 6 } ^b 
, 7 ) — >a 


^ I — 

f fit Y F3t Y fHt Y FF T Ffi 

fllgebralCalc Other PrgnIO|Clear a-z... 

: coe-f 
: Func 
: Loc^ 
: n.i 
: --^3^ 
006] 
: For 
: lefi 

;fCn) 
al i 

: -500 ,0,0,0,0,0j0>-J-a 

300, 10005, 16512, 19510,20020, 14515,6 

i, Ipfi 

:-(:augnent(: ^ "dotP^a, b)--"i> , a) , 7)->a 

I- coef f( 3000) 

3.9739E92eJ 

3.9739422655873E926 

MAIN RAD APPRDK FUNC 1/3:0 
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More Exciting Trig with the TI-92 

presented by Barbara Leitherer, 

4th International Symposium on DERIVE and the TI-92, Saeroe, Sweden 


I. Radian Measure - a FUN activity and concept 

In order to the "radian measure concept", the following toolbar menus need to be accessed from the geometry 
viewing window: 



The Curves and Polygons The Points and Lines Toolbar The Construction Toolbar 

Toolbar Menu to create circles Menu to construct a radius Menu to transfer measure- 

and arcs: i (segment): ments: ([ 



The Measurement Toolbar Menu to 
measure the length of an arc, a radius, 
a circumference, an angle and to do 
additional calculations: « 



The Display Toolbar Menu to make 
an arc graphically stand out from the 
rest of the circle: ¥ . 



1. Open the geometry window. 
On one side of the viewing 
window construct a circle. 
Then construct the radius of 
the circle (a segment) from 
any point on the circle to the 
center point of the circle. 


2. Measure the length of the ra- 
dius and the circumference 
of the circle and label each 
measurement properly: right 
after each numerical value 
appears, press A and the 
right cursor key > and type 
"= radius" and 

"=circumference". Then drag 
each label away from the 
circle area. 


3. Now access the construction 
toolbar menu and transfer 
the length of the radius to the 
circumference which creates 
a new point on the circle. 

Pay ATTENTION to the or- 
der in which this task has 
to be performed. Select first 
the endpoint of the radius lo- 
cated on the circle, second 
the circle, third the measure 
of the radius. 


4. Now overlay a minor arc from the new point constructed in step 3 and hook it up to the endpoint of the ra- 
dius which is located on the circle. Display the smaller arc in thick format, label it "=arclength" and drag text 
and numerical value away from the circle area. 
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5. Repeat the procedure from 
step 3 and prick off a second 
point which is a radiuslength 
away from the last point 
constructed. Continue to 
overlay the circumference with 
more radiuses. 

How may "radiuses" will ap- 
proximately fit around the total 
circle? 


6. Use the measurement toolbar 
menu and divide the 
circumference by the length of 
the radius. 

What number do you get? 

What does your answer mean? 

7. Now calculate the length of 
the remaining "little" arc. How 
many radians of the total cir- 
cumference does it represent? 


8. Now draw a radius to the point 
which was constructed in step 
3. How big is the central angle 
formed by the two "radii" of 
the circle? Give your answer 
in degrees, label it "=angle" 
and drag it to the right above 
the other measurements. 


9. Let's finish this project with a little exercise: 

A) Before you move on, fill the provided table with the measurements from the left and right of the circle. 



radius 

circum- 

ference 

arclength 

angle 

circumfer/ 

radius 

length of 
the little arc 

radians of 
the little arc 

original 

circle 








bigger 

circle 








smaller 

circle 









B) Drag the circle dynamically and make its radius once bigger, once smaller. Complete the missing 
information and answer the following question based upon your answers from the table: 


Which of the performed measurements DO, which DO NOT change in this exercise? 


II. Transferability to other topics 

1 . Defining Trigonometric Functions of any angle 


eg sin 0 


^ create a circle and measure its radius or use the measurement toolbar menu and select 
"equations and coordinates" 

^ create a point and its coordinates on the circle 
^ create and measure an angle 0 
^ find sin 0 

^ find ratio y/r (either through measurement toolbar or Data/Matrix Editor) 

^ change radius 

again observe sin 0 and y/r 



2 . 

3. 

4. 


Use for sign chart of different trigonometric functions in all four 
quadrants 

OR .... for Properties of sine and cosine: cos(-0 = cos {t) 

sin(-0 = - sin(t) 


SPACE for discussion 


II. References: 


1 . Classroom materials I used in my Preealeulus eourses taught at Carroll Community College 

2. TI-92 Guidebook, by Texas Instruments, 1995 

3. Preealeulus, Funetions and Graphs, by Demana, Waits, Clemens and Foley, 3rd edition, 1996. 
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WHAT ARE THESE PROBLEMS ABOUT? 

Alfonso J. Poblacion Saez 
E.U. Politecnica, Valladolid, Espana 

One of the main aspects I like to tackle in this section is to propose problems for which we 
can use DERIVE to solve (although I have not to be strictly necessary to do it) in order to 
reach different goals: to enjoy ourselves thinking and solving them, and to find the power of 
DERIVE programming and doing calculations, and why not, pointing out new ways to 
improve the system where it could appear more limited. So, I hope most of you (students 
included) may find useful or at least entertaining the problems proposed. 

This time I chose several questions with a common denominator which you also have to 
discover. Solutions and DERIVE files will appear in next newsletters. 

® What is the smallest four-digit number that belongs to the residue class 13 Zn and to the 
class 5 in Zu. And the biggest? 

(D In a livestock fair, a farmer bought a hundred animals which costed him 200 000 pts (the 
peseta is the Spanish monetary unit). Hens were valuated in 100 pts each, every pig in 
8 000 pts, and cows were sold at 40 000 pts each. How many animals of every sort did he 
buy? 

d) A mason had to floor three sqared factories with a certain number of tiles. In the first 
factory, he put the half of the total amount of tiles. For the second one he arranged a 50% 
more of rows than the third one. When finished, he found that he had one tile left. What 
was the total amount of tiles knowing that the factories were the smallest ones that verify 
all conditions referred? 

® A certain three-digit number yields a quotient of 26 when divided by the sum of its digits. 
If the digits are reversed, the quotient is 48. What is the smallest three-digit number for 
which this is possible? How many are three? 

® Could you prove using DERIVE that there exists a unique magic square (excluding 
symmetries) of order three that uses the first nine numbers? (A magic square is a square 
array of numbers in which the sum of each row, column, and main diagonals is identical). 


Fitting to thatACDC Contribution I can announce for the next DNL a TI-92 Program for Solving 
Diophantic Equations. It is a nice coincidence that it was also contributed by two Spanish DUG 
members, Leandro Tortosa and Javier Santa Cruz. See also Andreas Eder's contribution an page 33 
referring to J. Wiesenbauer's TITBITS 1 1 drom DNL#27. 


'' AC DC is Amazing (or Amusing Comer of the DERIVER' s Curiosities) is founded by Alfonso Poblacion. 
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150 Years of Pi's 250 Decimals 

Tomass Romano vskis 
University of Latvia, Riga, Latvia 



To demonstrate the power of computers and of DERIVE we can use the computation of n up to thousands of 
decimals up to 500000. The history of competitive efforts to find as many digits of n as possible began in the 
middle of the last century. In 1841 W.Rutherford published the number n with 208 digits in "Philosophical 
Transactions", 1841, Part 11, pag. 238. At the same time the Viennese mathematician S.Straschnizky asked the 
1 6 years old prodigy Zacharias Base to compute n using the formula: 


71 _1 1 

— = tan — h tan 
4 2 


1 

— + tan 
5 


\_ 

8 


To illustrate his capabilities Td like to mention a well known example of a multiplication. To find the product of 
79532853 * 93758479 took him 54 seconds. Within two months Base found the first 205 decimals of n and he 
sent his result to H.K. Schuhmacher, editor of the "Astronomischen Nachrichten". He compared Base's results 
with Rutherford's publication. The first 152 digits were equal, but from the 153th digit they did not even have a 
single digit in common. Rutherford had used another formula: 


71 , -1 1 -1 1 -1 1 

— = 4 tan tan h tan — 

4 5 70 99 


In that situation Schuhmacher wrote to Br. Clausen in Borpat (now Tartu in Estonia) and asked him to let him 
have his results, because it was known that he, Clausen, had calculated n for 250 places. He had even done that 
twice: once using the equation: 


71 = 8 tan ^ — + 4 tan ^ — 


and then using: 


71 = 1 6 tan ^ 


1 . -1 

4tan 

5 


1 


The values for tan ^ (1 / 3), tan ^ (1 / 5), tan ^ (1 / 7), tan ^ (1 / 239) and for n were published in the 
"Astronomi-schen Nachrichten" in 1847, Nr 589, S 207 - 210. So Schuhmacher found out that the boy prodigy 
had calculated n accurately to 200 decimals, and Rutherford's last 56 digits were wrong. 

Today we need not publish the values for arcustangens and for n again as everybody can now repeat the process 
using DERIVE. Set Precision Approximate and Bigits 250, edit Rutherford's, Base's or Clausen's formula, press 
one key and within fractions of a second you will get the 250 decimals places of Clausen. 

What now deserves more of our admiration? Rutherford, who made human mistakes, the Base boy, who made 
no mistakes, the human "super-calculator" Clausen or superquick perfect DERIVEl That is a matter of taste. 
Beep in my heart my favourites are T. Clausen and DERIVE. Using DERIVE you can go farther: by simplifying 
all formulae presented above in Exact Mode you will receive n exact in each case. Just try! If you approximate 
asking for 250 decimals the result will appear ten times faster. I assume that n in that case is computed by a 
different algorithm. 

DERIVE users - and not only they - should acquaint themselves with T. Clausen a bit more closely, because he 
was really a human super-calculator. He was born 1801 into a poor Banish family. At the age of 12 he could 
neither read nor write. A hobby astronomer, by name of G. Holst, gave him both some basic knowledge and the 
opportunity to attend school. When he was 20 he computed the planet orbits so well, that K. Schuhmacher and 
comet researcher V. Olbers took note of his talents. Later he published 150 scientific articles. He even solved 
problems which were too difficult for the great K. Gauss. But he became famous through his phantastic 
calculating abilities. He found errors in publications of world-renowned scientists and in logarithm tables 
prepared by 

K. Matissen. T. Clausen showed that the Fermat number 2^^ + 1 is not prime but can be factorized in 
274177 * 67280421310721. Check it immediately with DERIVE. 

Riga, August 27, 1997 


Tomass was delegate at the RJN-Conference at Saeroe in Summer 1997. We spent some time in 
talking and learning about Latvia. We all remember his lecture about Kepler's Clock and his 
presentation of a folk song of his home country. Tomass sent his artcle in German and I fear that my 
translation is not as excellent as his German paper. My best wishes to Riga, Josef 
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The New Supercharged TI-92 Plus Module with Advanced Mathematics Software 
for University Level Mathematics and Science 

by Bert K. Waits and Franklin Demana 

Department of Mathematics, The Ohio State University 

Version 12-26-97 

If you have suggestions please email them to: waitsb@math.ohio-state.edu 

Texas Instruments recently announced the TI-92 Plus Module, a powerful upgrade for the TI-92. This 
new module upgrade plugs into the back of any TI-92. The original TI-92 has built-in a versatile 
computer symbolic algebra system, a computer interactive dynamic geometry system, five graphing 
utilities including 3D graphing, matrix algebra, statistics, data, table, programming, and text editor 
functionality. The new TI-92 module upgrade contains considerably more memory (additional RAM 
of 128K AND a new User Data Archive memory space of 384K), a new Flash ROM that is 
upgradable, AND new university level mathematics, science, and engineering software. The 3D 
graphing feature has been improved by including contour plots. Implicit graphing of relations is now 
possible. Any implicit, contour, or 3D graph (in any style) can be rotated (“spun”) creating dynamic 
and exciting real-time animations and visualizations. Additional major new functionality includes 
numerical and exact analytic solutions of most 1st- and 2nd-order differential equations, exact solution 
of systems of linear and non-linear equations, advanced linear algebra, more statistics, a general 
numerical solver, units, and formatting improvements including removal of the circular definition 
error. Also assembly language programming is now possible. Some examples of the new “TI-92 Plus 
Module” functionality are illustrated in this paper. 

The new flash memory technology is a significant new feature. It makes the TI-92 electronically 
upgradable, adding considerable value to the product for students. With the TI-92 Plus Module, a 
student can download the latest software versions AND new applications direct from the TI web site 
(www.ti.com/calc). 

Section I contains examples and Section II contains a brief “workshop.” 

Section I Examples 

Differential Equations - Numerical Solutions 

Differential equations and initial value problems can be solved using the new module. 

Example I-l: A FIRST ORDER DIFFERENTIAL EQUATION 

Consider the standard logistic growth model given by the first order differential equation 
dP 

— = P^{C- P) where P =P(t) is a population at time t. Let k = 0.001 and C = 100. First we 

dt 

visualize the differential equation with a slope field (Figure 1). The slope field display occurs as the 
default view if no initial conditions are specified. 
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Figure 1. 


A slope field for 


dt 


[-10, 100] by [-10, 110] 
O.OOIP(IOO-P). 


If we specify P = 10 when t = 0, we obtain the solution to the corresponding initial value problem 
(Figure 2). 
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Figure 2. 


dP 

The solution to the initial value problem = 0.001P(100 - P), P(0) = 10 . 

dt 


Initial condition(s) can also be selected interactively using the menu F8 [IC], or by using lists of initial 
conditions in the editor. The table feature also provides useful numerical data as shown in 

Figure 3. 
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Figure 3. A graph and table of values for the so- 
lution to the initial value problem 


dt 


= 0.001P(100-P), P(0) = 10. 


Example 1-2: A SECOND ORDER DIFFERENTIAL EQUATION 
One standard model for the damped motion of a simple pendulum is 


6^{t) = -C6\t)-w^ ^mO{t)y^\iQVQ C is a constant representing friction and air resistance and 


w = 


. Here L is the length of the pendulum rod and 0(t) is the angular displacement of the 


pendulum rod from the vertical (equilibrium position). 













p 54 Bert Waits & Frank Demana: The SPAM D-N-L#29 


According to Dreyer this equation cannot be solve analytically (p 248, Modeling with ODE’s by T.P. 
Dreyer). However, it is an easy task for the new TI-92 Plus Module! We solve the problem for 
constants L = 2, C = 0.1, and g = 9.80665. The solution graph (Figure 4) is shown for the initial 
conditions yil= -3 and yi2 = 0. Note that we have turned on ‘‘labels” from the GRAPH FORMATS 
screen ([♦]] [F]). 



Figure 4. A solution orbit of the damped pendulum problem. 


The axes settings are set using the Style menu from the “Y=” Editor. Figure 5 shows several 
solutions using the interactive explore feature (Menu F8 IC from the graphics screen). One plot was 
set with yl = 0 and y2 = 5. The second was set with yl= 0 and y2 = 6 (just input these values after 
pressing the F8 key). 



Figure 5. Additional solution orbits of the damped pendulum problem using the new 

interactive initial condition feature. 

FOR FURTHER EXPLORATION 

For the damped simple pendulum, discuss the singularities (center and saddle) and their stability 
(stable or unstable). 

Differential Equations - Exact Symbolic Solutions 
Example 1-3 : A FIRST ORDER DIFFERENTIAL EQUATION 

We can use the new CSA (computer symbolic algebra) command deSolve( to find exact solutions to 
the most commonly studied 1st- and 2nd-order ordinary differential equations. The analytic solution 
to the logistic differential equation discussed earlier in Example I-l is shown in Figure 6. 










Figure 6. Using a new TI-92 Plus Module command deSolve to find an exact solution of a 
differential equation. Notice the general solution is found (@1 is the arbitrary constant) with ease. 



Figure 7. Determining the exact solution satisfy- 
ing the initial condition P = 10 when t = 0. 


In Figure 7 the regular TI-92 solve( command is used to find the exact solution satisfying the initial 
condition P = 10 when t = 0. Figure 8 shows the population graph and the value of the solution (ys, 
the population) at ^ = 40. 




[-10, 110] by [-10, 110] 

Figure 8. A graph of the exact solution and the exact solution when t = 40. 
















p 56 


Bert Waits & Frank Demana: The SPAM 


D-N-L#29 


Example 1-4: A CLASSIC SECOND ORDER SYSTEM OF DIFFERENTIAL EQUATIONS 

Predator-Prey problems are usually expressed as non-linear systems of coupled differential 
equations like those shown in the first panel of Figure 9. The solution graphs of the two interacting 
populations (as functions of time t) P1(T) = yl(t) and P2(t) = y2(t)) are also shown in the first panel of 
Figure 6. The window is [0, 10] by [0, 25], initial conditions yil= 2 and yi2 =5 when t = 0, and 0 < t 
< 10. The “thick” plot is yl (we did this using [F6: STYLE] from the Y= Editor.) Note: you should 
turn off the slope field option from GRAPH FORMATS screen. 



Figure 9. A graph of the exact solution and the exact solution when t = 40. 


Figure 9 shows the individual population graphs as function of time t. Now we change the AXES 
settings to obtain perhaps the more expected phase plane portraits (the orbits). We select “X Axis: yl 
and Y Axis: y2” from [F7], the AXES menu in Y^Editor application (choose 2: Custom as shown in 
the first panel of Figure 10). Figure 10 displays the results of graphing with initial conditions given by 
the two lists in yil and yi2 as shown in the second panel of Figure 9. 



Figure 10. A graph of the exact solution and the exact solution when t = 40. 


The second panel in the last row of Figure 10 shows a TRACE numerical result at time = 1.1 for one 
of the “orbits” (determined by what pair of initial conditions?). We have changes the ymax parameter 
to 50 in Figure 10. 

FOR FURTHER Exploration 

Approximate the center {a,b) of one of the orbits. Change the initial conditions to yl = a and yl = b 
and the F7 AXES settings back to “Time”. Before you replot can you predict what the two popiltion 
graphs (as functions of time t) will look like? Graph and then confirm what you see using analytic 
paper and pencil methods of calculus. 
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CONTINUED.... 

******************************************************************* 

This paper is 22 pages in length. You can download it (now in 2012!!) from 

http://www.math.osu.edu/-waits. l/papers/ti92spam.pdf 


/ am sure that this information about the SPAM is interesting for DERIVE-Users who haven't 
worked with a hand held tool until now. It is fascinating how many CAS-features can be 
implemented in such a machine. For me - and for many of us - the TI-92 and DERIVE are not 
separated of each other. In most cases it is possible to convert DERIVE activities into TI-92 
activities. And I often successfully try to transfer TI-92 ideas onto the PC using the 
advantages of working with more memory and a larger screen. Josef 


Additional comment (2012): 

You can download a couple of papers from the website given above. They are still of 
importance because TI-NspireCAS is heir and successor of TI-92 SPAM, TI-92 PLUS and 
Voyage 200. Bert Waits and Frank Demana were important propagators of CAS in the US - 
and world wide. I hope that Bert can read these lines. I met him at several occasions and am 
keeping wonderful memories in my mind, Josef. 


Keith Eames, UK keitheames@email.msn.com 

I have a slight problem with DERIVE for Windows, what must the settings be in order for circles to 
look like circles when plotted, similarly for coordinate geometry so that squares look like squares 
when printed. If this is in the Derive handbook then I apoligise but I left it at work! 

DNL: Answer from Soft Warehouse: 

Automatic setting of the aspect ratio so that circles always look like circles is not part of DERIVE for 
Windows version 4. Ox. We do plan address this issue in the next major release of DERIVE for 
Windows, version 5.0. In the meantime, you'll probably have the most success setting up horizontal 
and vertical grid values that work best for the size of the 2D plot window used. For instance on my 
screen, if you vertically tile one Algebra and one 2D-Plot window, then issue the Options Grids com- 
mand and enter a horizontal interval of 7 and a vertical interval of 8 a circle looks pretty circular. 
You can also manually resize the 2D-plot window using the resizing handles to further refine the dis- 
play of circles and squares. 

Yes, we have received similar complaints. In DfW 5.0, we plan to add the capability of saving 2D 

& 3D plot windows to TIFF, JPEG, GIF, etc. graphic output files. So, tell your friends to hang on! - 
Aloha, Theresa 

Gordon Evans 100302.2510@compuserve.com 

I apologise for this request, but the information I have on DERIVE is limited. How did you set about 
plotting fileds of vector functions such as: 

G(x,y) = - (i y + j x)/sqrt(x^2 + y^2) i and j being unit vectors. 


Can you help please? - Regards Gordon Evans 
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Neil Stahl, Menasha, USA nstahl@uwc.edu 

IVe used the following functions with DERIVE for DOS. I haven't tried them with DfW but probably 
you can use them to find a way to plot the fields there. 

It's a two-stage process. First GRID puts dots at the points, then VFIELD generates line segments for 
the arrows from those points. There are no heads on our arrows. In DfD one could set the parametric 
range for lots of curves at once. I hope that's possible in DfW. You should probably keep graph 
colors from cycling. 

E.g. you could type, then simplify and plot: 


GRID (-2, 2, 5, -2, 2, 5) 

VFIELD ( [-y,x] /sqrt (x^2+y^2) , -2 , 2 , 5 , -2 , 2 , 5 ) 

GRID generates a lattice of 2D points for a<=x<=b and c<=y<=d." 

GRID (a,b,m,c,d,n) : =VECTOR (VECTOR ([x,y],x,a,b, (b-a)/m),y,c,d, (d-c)/n) 

VFIELD can be used to plot elements of the 2D vector field v2 at points indicated as in GRID." 

VFIELD (v2, a,b,m, c, d, n) :=VECTOR (VECTOR ( [x, y] +t_^v2, x, a,b, 

(b-a) /m) , y, c, d, (d-c) /n) 

DFIELD plots the direction field u (probably depending on x and/or y) at points indicated as in 
GRID." 

DFIELD (u, a,b,m, c, d,n) : ^VECTOR (VECTOR ( [x, y] +t_^ [l,u] /SQRT (1+u^u) , 

x,a,b, (b-a) /m) , y, c, d, (d-c) /n) 

Terence Etchells and A1 Rich: 

We are experience a problem with DFW on our network. We are using 3.11 on pentium 
class chips with 16mb RAM. When DFW loads it uses 8mb of Ram and WIN32 uses just 
over 4mb!! (it does not do this in Windows 95). This makes the copying and pasting of 
graphics into, say Word, difficult with lack of memory messages in abundance. Is this is a 
problem of using WIN32 and is this problem likely to occur with any 32 bit prog in Windows 
3.11? As anybody else experienced this problem and is there a solution? We have tried all 
the suggestions on your FAQ on the web, it made no difference! 

Yes. The memory management is very inflexible when DERIVE for Windows (DfW) is running 
under Win32s which, in turn, is running under Windows 3.1x. 

Windows 3.1x is basically a 16-bit protected-mode MS-DOS program. Windows 95 is a far superior 
platform for running computer algebra systems, like DERIVE, because it is based on a true 32-bit 
operating system. 

Thus, I would highly recommend that DERIVE users seriously consider converting to Windows 95, or 
the forthcoming Windows 98. Note that version 5 of DfW will NOT be able to run on Windows 3. lx. 

- Aloha, Albert D. Rich 






